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ABSTRACT 

Massive, galaxy-scale outflows are known to be ubiquitous in major mergers of disk galaxies in the 
local universe. In this paper, we explore the multiphase structure and power sources of galactic winds 
in six ultr aluminous infrared galaxies (ULIRGs) at 2: < 0.06 using deep integral fleld spectroscopy with 
the Gemini Multi-Object Spectrograph (CMOS) on Gemini North. We probe the neutral, ionized, 
and dusty gas phases using Na I D, strong emission lines ([O I], Ha, and [N II]), and continuum 
colors, respectively. We separate outflow motions from those due to rotation and tidal perturbations, 
and flnd that all of the galaxies in our sample host high- velocity flows on kiloparsec scales. The 
properties of these outflows are consistent with multiphase (ionized, neutral, and dusty) collimated 
bipolar winds emerging along the minor axis of the nuclear disk to scales of 1 — 2 kpc. In two cases, 
these collimated winds take the form of bipolar superbubbles, identifled by clear kinematic signatures. 
Less collimated (but still high-velocity) flows are also present on scales up to 5 kpc in most systems. 
The three galaxies in our sample with obscured QSOs host higher velocity outflows than those in the 
three galaxies with no evidence for an AGN. The peak outflow velocity in each of the QSOs is in the 
range 1450 — 3350 km s~^, and the highest velocities (2000 — 3000 km s~^) are seen only in ionized 
gas. The outflow energy and momentum in the QSOs are difficult to produce from a starburst alone, 
but are consistent with the QSO contributing signiflcantly to the driving of the flow. Finally, when 
all gas phases are accounted for, the outflows are massive enough to provide negative feedback to star 
formation. 

Keywords: galaxies: evolution — galaxies: interactions — galaxies: ISM — galaxies: kinematics and 
dynamics — ISM: jets and outflows — galaxies: quasars 



1. INTRODUCTION 



et al. 



Massive, high- velocity outflows are present in almost 
all major galaxy mergers in the local universe ([ Heckman 
^'^^^ [2OOOI IRupke et al.||2002i |2005a|brct .Martin. 

Westmoquett e et al.||2012| [Soto Martin] 



2005 



2012 



1990 



2006 



Soto et al.|2012|) . These kpc-scale outflows exist in 



both the neutral and ionized gas phases. They have high 
velocities (typically a few hundred km s~^, though reach- 
ing >1000 km s~^ in some cases); large mass outflow 
rates (estimates suggest a few tens of percent of the star 
formation rate); and large extents (>5 kpc). They have 
been thought to be powered primarily by star formation 
dRupke et al.||2005at [Krug et al.||2010| ). Their properties 
are consistent with requirements lor negative feedback 
on star formation and black hole accretio n in mergers, as 



Springel et al. 



suggested by numeri cal simulations (e.g, 
^005] [Hopkins et al] [2005 ) . 

However, these conclusions rest on studies using long- 
slit optical spectroscopy of a single gas phase. Such stud- 
ies do not probe the detailed structures of winds, and 
have not compared in detail the properties of different 
gas phases. In particular, disentangling motions in the 
ionized gas (as probed by strong emission lines) due to 
rotation, tidal motions, and outffows is difficult without 
spatially-resolved data. The present ubiquity of sensitive 
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optical integral fleld spectrographs on large telescopes 
has made it possible to study the winds in these systems 
in much greater detail. 

Late-stage, major mergers of massive disk galaxies in 
the local universe are synonymo us with ultraluminou s 
infrared galaxies (ULIRGs; e.g., Veilleux et al. 2002). 
ULIRGs have been studied f or many years using integral 
fleld spectroscopy (IFS; e.g., Colina et al.||2005| and ref- 
erences therein). However, these early studies used 4m 
class telescopes at relatively low spectral resolution, and 
did not typ ically address the outflows in these systems 
(though see [Colina et al.||1999D . 

Recent deep observations on 8 — 10m class telescopes 
at high spatial resolution have made possible the de- 
tailed study of the multiphase wind structure an d power 
sources in ULIRGs (Arribas et al. 2008; Shih fc Rupke 



2010] [Rupke fc Veineuxf2CTl l ; West moquefte et al.|2012| . 
Shih Kupke[ ( [2010 ) conflrmed the outflow veloci- 
extent mea s ured previous ly in 



ties, mass, and radia 



F10565+2448 (Rupke et al. 2005c 



Martin 



2006), but 



also showed that the outflowing gas has a complicated 
structure that is inconsistent with lifting of r otating 
gas fr om the disk (Martin^ 2006} . Rupke & Veill eux] 
(2011) conflr med the detection o f an extended outflow 
in IVlrk 231 (Rupke et al. 2005a) and showed that this 
outflow emerges in all projected directions and is likely 
powered by the central QSO. This galaxy may host the 
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first known exan iple of QSQ-mode feedback in a rnajor 



merger. (See also Fischer et al. 2010a Feruglio et al.|201Q 



for molecular gas detections of the outflow. ) Westmo- 
quett e et "aT] (2012) presented the first large-scale 11^'IS sur- 
vey ot outflows in ULIRGs. These authors focused on the 
ionized gas phase, and detected ubiquitous kiloparsec- 
scale outflows through the presence of broad, blueshifted 
velocity components. The high velocities they detected 
in some systems (reaching several thousand km s~^) im- 
plicate an AGN in powering the outflows. 

In this paper, we present an IFS study of six ULIRGs. 
This study differs from early ULIRG studies in its high 
depth, spatial resolution, and spe ctral resolution. It 
also differs from the recent study of |Westmoquette et al. 
( [2012| ) by our ability to probe multiple gas phases (neu- 
tral and ionized) at higher spatial resolution (enabled 
by 0^2 spatial sampling) and by our modeling of the 
resulting wind structure. The data on two of these 



systems were first presented in Shih & Rupke (2010) 



and Rupke & Veilleux ( 201l] ), and we re-reduce and re" 
analyze these data along witn four new systems. Five of 
these systems were first su rveyed for ou tflows using long- 
slit spectroscopy (Heck man et al.|2000 Rupke et al.^2002] 
|2005a|c| [Martin 2005). 

In Section [2j we discuss our observations and data 
analysis. We discuss our analysis procedures and wind 
modeling in SectionjS) We then detail the properties of 
individual objects in Section|4] The results on individ- 
ual systems are synthesized to form a sample-wide pic- 
ture in SectionjSj with the crucial results found in Fig- 
ures |49l - [56l We summarize in Section[6l 

2. SAMPLE, OBSERVATIONS, AND REDUCTION 

The six systems in our sample (Table [l]) were se- 
lected for their proximity and high infrared luminosity: 
z < 0.06 and log(LiR/L0) > 12. All six are members 



of the Rev ised Bright Galaxy Sample (RBGS; [Sanders 
et al.||20Q3 ), which was assembled based on Infrared As- 
tronomical Satellite data. They span a range of nuclear 
properties and merger classes, though all are thought to 
be in the l ate stages of a major merger (mass ratio < 3; 
V eilleux et a l. 2002). 

I'he galaxies in our study were observed with the Inte- 
gral Fiel d Unit in the Gemini Mult i- Object Spectrograph 



(GMOS; |Allington-Smith et al.] |2QQ2 
on the Gemini Worth telescope. In 



eac 



Hook et al.||2004 ) 

1 case we used the 
integral field unit in 2-slit mode, with the B600 grating 
and r' filter, yielding wavelength coverage from 5600 to 
7000 A and a spectral resolution averaging 1.6 A at the 
central wavelength of 6300 A. (This spectral resolution is 
measured directly from night sky emission lines in each 
spaxel.) 

To assess image quality and relative photometry in 
galaxies with an identifiable nucleus, we fit a broad Gaus- 
sian to the nucleus in each exposure. Variations in im- 
age quality and cloud cover between exposures led to at 
most a 30% difference in this peak flux between the high- 
est and lowest fluxes. A small number of exposures with 
larger variations were discarded. Delivered image quality 
through the IFS was difficult to assess for most sources, 
since there are no point sources in the field of view of 
the integral field unit (except for the QSO in Mrk 231, 
which has a full width at half maximum, or FWHM, of 
0^5 — 0^6). However, fits to multiple field stars in 12 



of the 13 acquisition images in our dataset yield a mean 
FWHM of 0V59, with a standard deviation of 0^16 and 
a range of 0V4 - 0^9. 

The native FOV of GMOS in 2-slit mode is 5" x 7". 
All sources were observed in two pointings. In four sys- 
tems, GMOS was dithered along the short axis, with 
some overlap between the two pointings. In the other 
two (F08572+3915 and VV 705), it was dithered between 
two nuclei. 

Table [2] lists further details of the observations of each 
source. 

We reduced the data using IRAF vl.l4 and vl.lO of the 
data reduction package provided by the Gemini Obser- 
vatory. We modified the Gemini software slightly to ac- 
commodate error propagation, and supplemented it with 
custom IDL routines. These IDL routines perform the 
following steps: (1) correcting single- wavelength images 
to a common center (due to differential atmospheric re- 
fraction); (2) applying a small shift to the wavelength 
calibration using the line centers of atmospheric emis- 
sion lines ([O I] 6300, 6364 A or Na I D 5890, 5896 A); 
(3) telluric correction (using a scaled version of the high 
SNR Mrk 231 nuclear spectrum as a template); (4) spa- 
tial rebinning of spaxels; and (5) mosaicing of dithered 
positions. 

The native 0V2-diameter hexagonal spaxels of GMOS 
are resampled to 0^1 square spaxels by the GMOS IRAF 
pipeline. To Nyquist sample the seeing disk while max- 
imizing sensitivity, we binned the O'.'l square pixels into 
0^3 square pixels. 

The relative flux calibration of our data is robust. 
However, flux calibration uncertainties at the red fllter 
transmission edge may affect the [N II] /Ha ratio at the 
<50% level in F08572+3915, which has the highest red- 
shift in our sample and thus an [N II] line near the fllter 
edge. 

The flnal result is eight data cubes (one per system 
dithered along the GMOS short axis, and two per system 
dithered between nuclei). The resulting cubes contain 
hundreds of spectra. 

3. ANALYSIS AND MODELING 
3.1. Identification of Nuclei 

The nuclei in these systems are heavily obscured, and 
the optical continuum peak is often not coincident with 
the underlying stellar mass peak. The stellar mass peak 
can be approximated by the location of the near- infrared 
(NIR) luminosity peak (e.g., jScoville et al.|[2QQQt . 

Archival HST F814W and F16UW images exist for each 
galaxy in our sample. To identify the nuclei within the 
GMOS data cubes, we located the NIR peak or peaks 
in each HST F160W image. We registered the NIR im- 
age with the F814W image by matching the positions of 
1 — 3 relatively unobs cured nuclear star clusters, as in 
Iwasawa et al. ( |2011b ). We then convolved each F814W 
image with a 5-pixel Gaussian to match the average see- 
ing in our data (0^6). Finally, we matched the convolved 
F814W and GMOS continuum peaks. The GMOS con- 
tinuum is slightly bluer than the F814W continuum {R 
vs. / band), but we do not expect the continuum peak 
to shift signiflcantly over this wavelength range. The re- 
sulting nuclear locations in the GMOS data are accurate 
to ^OVl, or a third of a spaxel. 
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The exception to this procedure is Mrk 231; the nucleus 
for this galaxy is easily identifiable as the bright QSO. 

3.2. Continuum, Emission Line, and Absorption Line 

Fitting 

We modeled the continuum and emission lines in each 
spectrum using UHSPECFIT, a suite of IDL software 
which optimizes continuum and emission li ne fits to spec- 
tra. UHSPECFIT is described elsewhere (Rupke et al 



2010). In brief, it masks emission line regions, tits the 



contmuum, and then simultaneously fits all emission lines 
in the continuum-subtracted spectrum. The fit is im- 
proved by repeating this process based on the parameters 
derived from the first fit. 

We fit simple polynomial continua (of order 3 — 4) to 
each spectrum. For Mrk 231, we also fit a scaled version 
of the extremely bright nuclear QSO spectrum ( [Rupke 
|fc Veilleux|2011p . The QSO point spread function (PSF) 
bleeds into many spaxels and must be removed. 

The continua in our galaxies are dominated by young 
stars, which have weak absorption lines in the wavelength 
range of our data (>5200 A rest frame). The stellar 
absorption line that would most affect our emission line 
fits is Ha; because our spectra do not cover other Balmer 
lines, we cannot constrain the strength of Ha absorption 
using other lines. However, observations of ULIRGs show 
that the equivalent width of H/3 absorption is typically in 
the range 3 — 7 A within 2 kpc of the nu clei of ULIRGs, 
and only 1 or 2 A higher at larger radii (Soto & Martin 



20101). 



The Ha emission line equivalent widths in our data 
cover a broad range (from —1 to —300 A), but typically 
dip below —10 A only at radii >2 kpc, and sometimes not 
at ah. Some systems (Mrk 273 and VV 705:NW) have 
Ha emission line equivalent widths of < —20 A in virtu- 
ally every spaxel. The equivalent widths vary spatially 
in a complex way, but generally increase with increas- 
ing radius. Furthermore, the Balmer absorption lines 
in ULIRG s are typically much b roader than the emis- 
sion lines (Soto & Martin 2010), minimizing their im- 
pact on the emission line tittmg. The most likely effect 
of the lack of Balmer absorption correction in our data 
would be an overestimate of the [N II] 6583 A/Ha and 
[O I] 6300 A/Ha emission line ratios 2 — 3 kpc from the 
nucleus in a few systems. Conservatively, this overesti- 
mate is no worse than 50% in those spaxels with emission 
line equivalent widths near —10 A. Given the relative 
widths of the absorption and emission lines, the likely 
effect is smaller tha n this. This error does not affect the 
conclusions of § |5.2[ 

The emission lines visible in the wavelength range 
of our data are Ha, [N II] 6548, 6583 A, and 
[O I] 6300, 6364 A. We fit 1 - 2 Gaussian velocity com- 
ponents to these emission lines in each spectrum in all 
cases but Mrk 273, for which we allowed up to 3 compo- 
nents. We fitted a common velocity and FWHM to all 
emission lines in a given velocity component. The broad 
emission line component in F08572+3915:N W req uired 
a more constrained fit, which is described in § 4.1.2| The 
required number of Gaussian components in each spaxel 
was determined by a combination of software automation 
and manual inspection. UHSPECFIT discarded all com- 



ponents that were narrower than the spectral resolution, 
too broad to be robustly distinguished from the contin- 
uum, or feh below SNR = 3 in the Ha and [N II] 6583 A 
lines. The exception is Mrk 231, for which we required 
SNR > 7 in Ha or [N II] 6583 A. 

For spaxels with multiple velocity components, we or- 
ganize them principally by FWHM into a (1) narrow 
and (2) broad component (except for Mrk 273; see §4.4). 
However, we adjust the classification of some spaxels by 
hand to achieve continuity in the velocity maps, partic- 
ularly in cases where the FWHM of the two components 
are similar. 

For five of the six galaxies in our sample, we found 
a >3cr component in Ha or [N II] 6583 A in >90% of 
spaxels. In the sixth galaxy, F08572+3915, we detected 
emission in >60% of spaxels. 

We fit the Na I D a bsorption lines using the method of 
Rupke et al. (2005b). This method robustly separates 
optical depth and covering factor, even for a blended 
doublet and multiple velocity components, by using a 
Gaussian in optical depth. The resulting values for opti- 
cal depth vary from 0.1 to 5 (the latter is the upper limit 
of values to which we are sensitive), but are typically of 
order unity. The covering factor can vary between and 
1, but is typically 0.5. These results are consi stent with 
previous studies of nuclear sp ectra of ULIRGs (|Heckman 
erar][2QQQ{ [Rupke et al.][2QQ2[ I2005b| . 

One velocity component m IN a I D provided good fits 
to spectra throughout the sample. For most systems, we 
ignored the nearby He I 5876 A emission line because of 
its weakness. The exceptions are VV 705 and Mrk 231; 
in these cases we fixed the line width and centroid of He I 
in each spaxel to the values determined by fits to other 
emission lines in that spaxel. 

The stellar contribution to the N a I D absorption lin e 
in ULIRGs is generally quite small ( Rupke et al.|2005'b ). 
In previous studies, the stellar contribution to IN a I D has 
been estimated by scaling the strength of Mg I b ( Heck- 
man et al.||2000[ jRupke et al.| [2002, 2005b). The current 
data do not probe Mg I b, or indeed other strong stellar 
features, so we cannot estimate this contribution on a 
spaxel- by-spaxel basis. For four of the six galaxies in our 
sample, measurements of nuclear spectra fr om previous 
work (Kim et al. 1995; Veil leux et al.L99 9a^ or from the 
SDSS allow us to estimate this contribution for at least 
one nucleus (F08572+3915:NW, F10565+2448, Mrk 273, 
and VV 705:NW). We find that the equivalent width of 
Mg I b is 1.0 A or less, yielding a stellar Na I D equiva- 
lent width of 0.5 A or less. Compared to the measured 
total Na I D equivalent widths in these sources, then, the 
stellar Na I D represents only 5 — 10 % of the line. 

3.3. Maps and Velocity Distributions 

The immediate results of our emission- and absorption- 
line fits to the data are maps of physical properties in 
each system. The first set of maps for each system (e.g.. 
Figure [T| displays continuum properties: (1) two-color 
Hubble^pace Telescope images (using Advanced Cam- 
era for Surveys Wide Field Camera F814W and 
data in most cases); (2) continuum fiux ratio (color) im- 
ages using HST data; (3) continuum images from our 
CMOS data; and (4) the rest-frame equivalent width of 
Na I D in each spaxel. 
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The second set of maps (e.g., Figure |3| shows the hne- 
of-sight (LOS) velocity of each emission- and absorption- 
hne component in each spaxel. We display the centroid 
velocity, the most blueshifted velocity (for the broad 
components), and the FWHM (with the spectral reso- 
lution subtracted in quadrature) of each component in 
each spaxel. Our line profile modeling is based on Gaus- 
sian velocity profiles, so the most blueshifted velocity is 
defined by the properties of the normal distribution: 

v^Q% = center of Gaussian profile (1) 
(50% of gas has lower outflow velocity) 

'^98% ^50% - 2cr (2) 

(98% of gas has lower outflow velocity) 

These definitions assume that negative velocities are 
blueshifted, and apply to individual Gaussian compo- 
nents in cases where we fit more than one velocity com- 
ponent. 

The sole exception to these definitions is Mrk 273. Be- 
cause Mrk 273 shows significantly redshifted emission as 
well as blueshifted emission, we define Vgg% for redshifted 
components in this system by adding, instead of sub- 
tracting, (J. 

The third set of maps (e.g.. Figure |4| presents Ha flux 
and two emission line flux ratios ([N II] 6583 A/Ha and 
[O I] 6300 A/Ha) of each component in each spaxel. 

The statistical properties of the derived LOS velocity 
fields are presented in two ways. First, we show his- 
tograms of the distributions of all LOS velocities in each 
system (e.g.. Figure [s]). Second, we tabulate the mean 
and extreme values of FWHM, v^q%^ and Vq8% in each 
galaxy with an outflow in Table |3[ In the latter, we in- 
clude only outflowing components in each spaxel. 

3.4. Structural Modeling 

A primary goal of this study is to study the struc- 
tural properties of the outflows in individual systems. 
Structural models enable us to better quantify the ex- 
tent, mass, momentum, and energy in an outflow, and 
thereby its impact on the galaxy. 

Previous surveys of winds in ULIRGs used a spheri- 
cally symmetric, mass conserving free wind to estimate 
mass, momentum, and energy in outflows | Heckman 
et al.||2000[ jRupke et al.||2002[ |2005"cl |Soto et al. 2012). 
I'hese studies made assumptions about the radial extent, 
radial distribution, and variation of velocity with posi- 
tion of outflows in these s ystems. Recent IFS studies 
of outflows in two ULIRGs ( |Shih fc Rupke||2010| |Rupke 
fc Veilleux||2QTT ) were able to mitigate these uncertam- 
ties by providing lower limits on the radial extent of the 
gas and detailed maps of the gas mass and velocity as a 
function of position. 

In this study, we employ two strategies to model the 
structure of the outflows in each system. First, we esti- 
mate mass, momentum, energy, and their flow rates using 
a single radius, mass conserving free wind. Next, we use 
the detailed IFS data to motivate alternative geometries. 
In particular, we use superbubble models to fit the data 
for two systems (Mrk 273 and VV 705) that show strong 
evidence for such structure. Other geometries are possi- 
ble (e.g., biconical outfiows, filled-in geometries instead 
of shells), but require modeling complexity beyond the 



scope of the current paper. 

3.4.1. Single Radius Free Wind (SRFW) 

For the Single Radius Free Wind (SRFW) model 
we us e the ti me- averaged thin shell formulae of |Rupke 



let al.| (|2005c|), which are derived by assuming a mass- 



conserving how. These formulae were modified by |Shih 



Rupke (2010) for IFS data by letting each spaxel (cor- 



responding to a unique angular coordinate with respect 
to the wind's center) have its own mass and velocity. 
These formulae are equivalent to computing the mass, 
momentum, and energy at each angular location on the 
thin shell, and then computing outfiow rates by dividing 
by the dynamical time tdyn = Rwind/"^ of each spaxel 
(or angular coordinate). We assume a single radius for 
the wind in each system, estimated from the maximum 
extent of the wind in the plane of the sky. Within each 
system, we use the same wind radius for the ionized and 
neutral gas phases. We use the velocity measured in each 
spaxel, deprojected using the assumed radius, as the true 
wind velocity at that angul ar location i n the w ind. 

Most of the fo rmulae in Rupke et al. (2005c") and' Shih] 
& Rupke (2010| are expressed in terms or column density, 
and are in this form only applicable to the neutral gas 
in our study. However, by computing the ionized gas 
mass in each spaxel from the Ha luminosity, and dividing 
by the dynamical time for each spaxel, we can compute 
momentum, energy, and their outfiow rates in each spaxel 
from their dependence on mass and velocity. 

The neutral gas mass in each spaxel is proportional 
to the H column density and the deprojected solid an- 
gle at that location (i.e., the solid angle subtended by 
the spaxel with respect to the outfiow center). To com- 
pute the H column density, we first find the Na column 
density, which depends on the measured optical depth 
and velocity width in a simple way for our assumption 
of a Gaussian velocity distribution in eac h component 
(see, e.g., the formula in [Rupke et al.|2002 ). This is con- 
verted into H column density using typical Milky Way 
dust depletion for Na (—0.95 dex; [ Savage fc Sembach 
1996|) a nd metal abundance of Na (—5.69 dex; 'lavage 

Sembach 1996), and an uncertain ionization fraction 
[y = 1- 7V(Na l)/A/'(Na) = 0.9; [Rupke et al.||20Q5bl ). 
A^(H) depends on y as (1 — y)~^ . The solid angle sub- 
tended by each spaxel is easily computed in our spherical 
model, and we then multiply this solid angle by the mea- 
sured covering factor of the neutral gas, to account for 
dumpiness in the outfiow. (The covering factor in each 
spaxel equals the fractional area of the spaxel that is 
covered by absorbi ng gas, and is a para meter in the ab- 
sorption line fit; see Rupke et al.|2005b for more details.) 

The ionized gas niass in each spaxel is computed from 
the Ha luminosity, assuming Case B recombination and 
T = 10^ K. It depends on electron density as n~^. We 
assume Ue = 10 cm~^, which is consistent with limits on 
np, for the superbubble in NGC 3079 ( n^, = 5-100 cm"^; 
Veilleux et al.|[ 1994} [Cecil et al.|2001| ). 

'I'he assumed radii in the SKl^ W model for each galaxy 
with an outfiow are listed in Table [U This table also lists 
the deprojected covering factor over the surface of the 
outfiow, as seen from the galaxy nucleus, assuming the 
given radius, and summed over all spaxels with outfiow- 
ing gas. This integrated covering factor (which we label 
^I/Att) represents the fraction of sightlines, for an ob- 
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server at the galaxy center, that wih encounter the wind. 
It differs from the opening angle of the wind, which de- 
scribes the solid angle contained by the outer boundary 
of the wind. The integrated covering factor is a lower 
limit to the wind opening angle. 

3.4.2. Bipolar Superbubble (BSB) 

The ionized gas kinematics in Mrk 273 
neutral gas kinematics in VV 705 :NW ( 



4.4.3) and 



4.5.2p suggest 



the presence of a Bipolar SuperBubble (BSBJTTollowing 
similar treatments of other nuclear superbubbles (e.g., 
Veilleux et al.'1994 ) , we used two prolate ellipsoidal bub- 
bles on either side of the galaxy nucleus. Each bubble 
has semi-principal axis lengths Rz (along the axis shared 
by the two bubbles) and = Ry = Rxy (perpendicu- 
lar to the z axis). We assumed that the gas lies on the 
bubbles' surfaces and that each parcel of gas is moving 
radially with velocity v = Vz{r/rzY^ where r is the dis- 
tance from the center of the wind, and Vz and Vz = 2Rz 
are the velocity and distance along the z axis. 

For Mrk 273, we optimized seven of the eight BSB pa- 
rameters {vz^ Rxyi Ftzi liiiG of sight inclination, position 
angle, and center coordinates) by minimizing the RMS 
difference between the model and the data. We found 
that it was not possible to robustly constrain the power 
law index n beyond the statement that n > 1 is favored. 
The best-fit z-axis velocity Vz varied with n, such that 
larger values of n favored larger Vz- However, the best 
fit values of other parameters varied little or not at all 
with n. 

For two different plausible values of n based on studies 
of other sources ( 2 and 3; see, e.g., 'Veilleux et al. 1994 



Cecil et al.||^2QQl| , we list in Table [5] the best tit B'St 



parameters tor the ionized gas in Mrk 273. In the best 
fit models, 67 Nyquist-sampled spaxels lie completely in- 
side the bubble edges. However, each spaxel has two 
outflowing kinematic components (the near and far sides 
of the bubble), leading to 134 data points in the fit. The 
RMS difference between the model and data for the red- 
shifted and blueshifted components (with the redshifted 
and blueshifted RMS values combined in quadrature) is 
~200 km s~^. We estimated the errors in each parameter 
by fixing all other parameters and varying the parameter 
until the RMS increased by 50 km s~^. 

For VV 705 :NW, we performed a similar procedure. 
Because this disk is almost face-on (§4.5.2), any bubble 
emerging along the minor axis of the disk will also be face 
on. The outflow in this system was detected primarily 
in absorption, so we could not detect the far side of the 
bubble (behind the disk). The resulting 38 distinct data 
points (in the best fit model) made it infeasible to fit all 8 
parameters of the BSB model. Thus, we fixed the cente r 
of the outflow to be the NIR luminosity peak (§4.5.2). 
We also constrained the cylindrical radius of the outttow 
in the plane perpendicular to the z-axis to be Rxy = 
1.0 kpc after inspecting the extent of the outflow in the 
plane of the sky. Furthermore, as with Mrk 273 we could 
only determine that n > 1. Thus, for n = 2 and n = 3 
we fit four of the eight BSB parameters {vz, Rz^ line of 
sight inclination, position angle). 

To compute the neutral gas mass in the BSB model, we 
measured the column density in each spaxel and multi- 
plied by the projected spaxel area. The ionized gas mass 
was computed as in the SRFW model, and momentum 



and energy followed from the deprojected velocity of each 
spaxel. Outflow rates come from dividing by the dynam- 
ical time, 2Rz/vz- 

3.4.3. Results 

The mass, momenta, and energies and their outflow 
rates derived from the SRFW and BSB mode ls ar e listed 
for each system and each gas phase in Table 6^ 



Some uncertainties remain in these models. First, 
there is still LOS uncertainty in the gas position. Sec- 
ond, we employ only "thin shell" models: we assume that 
the gas is confined to a small radial range {AR/R << 1). 
Thick shell models are like ly to yield smaller m asses, mo- 
menta, and energies (e.g., fRupke et al ||2005cl). Finally, 
there is uncertainty due to the ionization correction re- 
quired to co nvert the amount o f Na I t o the total amount 
of Na (e.g., [Rupke et al.||2002[ |2005c| . We assume that 
the Na I repres ents 107o ol the tot al Na, based on Galac- 

Since Na resides pri- 



tic sight lines (Rupke et al. 



marily in gas well-shielded 



2002) 



Tom UV radiation (because 



of its low ionization potential), the correction is likely 
not much larger than this. 

4. INDIVIDUAL SYSTEMS 

In this section, we discuss the power source and kine- 
matics of each system in detail. 

4.1. F08572^3915 

This system is a close pair that is likely near second 
pericenter. The system's energetics are dominated by a 
buried AGN in the NW nucleus, based on copious ev- 
idence. We review this evidence and then discuss the 
outflow in this system, which shows the highest velocities 
in our sample. Figures [l]-[9] present maps and velocity 
distributions for each nucleus. 

4.1.1. Power Source 

The NW nucleus is the source of the prodigious far- 
infrared (FIR) luminosity, as it completely dominates the 
radio and mid-infrare d (MIR) luminosit y of the system 
(Condon et al. 1990l j Soifer et al.]|2000D . Inspection of 
Spitzer Space Telescope archival images at 3.6/im and 
8.0/im shows that the SE nucleus contributes less than 
1% of the system luminosity at these wavelengths. The 
dominance of the NW source is consistent with the ~2 x 
10^ McT) of molecu lar gas detected in the NW nucleus 
( Evans et al.|2002 ); only upper limits are obtained in the 

Clear evidence for AGN in this system is absent in 
ground-based op tical spectra. Previo us optical sp ectral 
types of LIN ER ( [Veilleux et al.|1999aD or Seyfert 2 ( [Yuan 
et al.|[2QTQ| ) in each nucleus were based on shallow spec- 
tra that did not s ignificantly detect [O III] 5007 A or H/3 



dKim et al.||T998[ ). The SDSS spectrum of the NW nu- 
cleus is much deeper than these previous spectra, with a 
similar spatial footprint. An analysis of the SDSS spec- 
trum using our stellar continuum and emission line fitting 



^ We note that the mass, momentum, and energy outflo w rates 
for Mrk 231 are a factor of ~2 lower than those hsted in jRupke] 
|fc Veilleuxl(|2011| >. This is because in [Rupke V eilleux ; 2011| > we 
did not multiply the neutral gas quantities m each spaxel by the 
covering factor of the absorption line gas in that spaxel. In this 
paper we do include this factor. 
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software (§[2| yields Ay = 2.14 . and a spectral typ e that 
is borderline between H and C ( Kewley et al.[2QQ6 ). This 
spectral type is not truly nuclear however, since the line 
ratios vary significantly in the E-W direction across the 
SDSS aperture. Some of these extended line ratios are 
consistent with AGN-excited gas, though not uniquely so 
in the absence of [O III] data (Fig. [4]). The optical ap- 
pearance in HS T images also shows Complex extinction 
and scattering" ( [Surace et al.||1998| Fig. m). In contrast 
with the NW nucleus, the nucleus is clearly of spec- 
tral type H (Fig. [s]), even in the absence of [O III] or H/3 
data, consist ent with its optical morphology (a "typical 
galaxy core" ; [Surace et al.|1998 ). 

In the NW nucleus, IN IK recombination lines point to 
heavier obscuration th an in the optical: Ay = 6.7, based 
on Paa and Br7 fiuxes (Veilleux et al.|1999b|) and assum- 



ing Case B, Ue 10^ cm" 2, and T 10^ K (Osterbrock 
|l989). This nucleus is unresolved in the NIR with HSl] 
and it becomes more prominent with respect to the ex- 
tended host between 1.1 and 2.2/ im (the half-light ra dius 
decreases from 830 pc to 140 pc; Scoville et al.|20QQ). L- 



band and MIR spectra point to a very high Ay, based on 
deep absorpti on featu res from carbonaceou s and silicate 
dust (Dudley fc Wynn- Williams 1997; Imanishi & Dud 



2QQ2,.Armus et al ,20071). The MIR 



|rey|2000^. 

source radius is <13 peat 12.5jtim (|!SQiter et al. 2000) 
and <60 pc at 18/im (Imanis hi et al.|201 l'). The increas- 
ing inferred visual extinction and decreasing source size 
as tracers move to longer wavelengths point to a highly 
obscured, compact source. 

The most likely candidate to power the NW source 
is a heavily obscured AGN, as confirmed by MIR spec- 
troscopy. The infrared continuum shape and lack of PAH 
emission provide evidence that a deeply buri ed AGN 
powers ^70 % of the total infrared luminosity (Veilleux 
et al. 2QQ 9 |. Hot d ust emission also points to an AGN 



( [Armus et al. 2007). This source is barely detected by 



the Chandra X-ray Observatory^ though its hardness ra- 
tio is in principle consistent with heavy obscuration of 
an AGN ( ,Teng et al. 2009) . 

4.1.2. Kinematics 

The SE galaxy reveals modest non-circular motions, 
indicative of tidal dynamics, superimposed on ordinary 
rotation of an inclined disk. The line of nodes is uncer- 
tain, but the major axis of the bulge is 130° E of N. The 
projected peak-to-peak velocity amplitude along this line 
is ~150 km s~^. The appearance of the spiral arms in 
the HST image and the redshifted emission in the SE 
point to the near side of the disk lying in the NE. No 
absorption-line gas or emission-line outfiow is evident in 
this galaxy. 

The narrow ionized gas component in the NW galaxy 
is undergoing orderly rotation, again with modest non- 
circular motions mixed in. The central isovelocity con- 
tours are not perpendicular to the line of nodes, which 
we estimate to lie at 120° E of N if it is to pass be- 
tween the projected rotational maxima (±120 km s~^). 
If the maximum possible disk velocity were 150 km s~^ 
(200 km s~^), then the disk inclination would be 53° 
(37°). The midpoint of the rotational maxima is offset 
~0.5''NE of the nucleus. If the tidal features near the 
nucleus are the elongations of spiral arms, then the near 



side of the galaxy is to the SW of the nucleus. 

The broad emission line component in the NW nucleus 
of F08572+3915 required special treatment. We were un- 
able to independently constrain the line fiuxes and width 
of the broad, outfiowing gas component in each spaxel be- 
cause of the large width of this component, which caused 
severe blending between Ha and [N II]. The outfiowing 
component of the isolated line [O I] was too weak to pro- 
vide an independent constraint on the line width. We 
instead used the single aperture Sloan Digital Sky Sur- 
vey (SDSS) spectrum of the NW nucleus to disentangle 
the line widths and fiuxes of the outfiowing component. 
While our spectrum does not have the wavelength range 
to include [S II] 6717, 6731 A, the SDSS spectrum does, 
and the outfiowing component appears in this line. In 
both the narrow and broad emission components of the 
SDSS spectrum, we fit simultaneously all emission lines 
using a common galactocentric velocity and line width 
in each component. We constrained the [S II] line ratio 
to lie within its expected values based on atomic physics. 
The Ha, [N II], and [S II] lines together provided enough 
leverage to simultaneously fix the line width and line 
fiuxes of the outfiowing component. We then used these 
values to fix the [N II] /Ha ratio and line width across the 
outfiow in our data, and we fit only the central velocity 
and absolute fiux in the broad outfiowing component in 
each spaxel. Fig. [lO] shows example fits to the broad line 
in our data. 

The broad ionized gas component reveals the high- 
est velocity outfiow in our sample, with ('^98%) = 
—2800 km s~^ and a maximum value for Vgg% of - 
3350 km s~^. It peaks in velocity directly between the 
peaks in projected rotational velocity, and shows lower 
velocities 2 kpc to the S/SW (though still in excess of 
'^98% = —2000 km s~^). The blueshifted component is 
very broad (FWHM = 1500 km s~^), but it is not an 
AGN broad line region because it appears in other for- 
bidden lines ([S II]) and is spatially resolved. 

The neutral outflow in the NW galaxy is located 1 — 
2 kpc SE of the nucleus, and reaches velocities up to 
-1150 km s-^ (with (^93%) = -900 km s-^). 

The outflow in F08572+3915 is not well enough re- 
solved for a structural model, but the data offer some 
clues. The emission-line flow extends along the galaxy 
minor axis, though in a direction behind the near side 
of the disk. However, at this point in the disk the red- 
dening is low. The decrease of outflow velocity in this 
direction is consistent with the far-side of a minor-axis 
outflow. The other direction along the minor axis shows 
heavy reddening; the outflow is perhaps obscured along 
this direction. The obscuration along this direction may 
also indicate a dusty flow. However, no neutral gas is de- 
tected here; the neutral gas outflow instead appears SE 
of the nucleus along the line of nodes. The neutral gas 
partly overlaps the emission-line outflow, but its overall 
morphology is unresolved, and it is unclear how the two 
are related. 

4.2. Fl 0565^ 2448 
The IFS data for this galaxy, a triple system dominated 



in luminos i ty by one member, was flrst presented in Shih 
& Rupke (2010). We review the nature of the power 



source m this system and revisit the interpretation of the 
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IFS data based on a new reduction and analysis with the 
current data pipehne and a finer binning (0^3 spaxels, vs. 



ff!6 in iShih & Rupke||2QlQ[ ). Figures [n]-{T5] show m 
elc 



and velocity distributions tor this system. 



4.2.1. Power Source 



Most of the luminosity in this triple, at all wavelengths, 
arises in the westernmost nucleus. A second and third 
nucl eus are located 8^.^Q SE ([Scoville et al.|2QQQD and 25^8 
E ( iMurphy et al.| 1996') of the W nucleus 'i'he E and SE 



nuclei are not detected in 1.49 GHz VLA images d Con- 
don et al.|1990 ) or Chandra X-ray images (Iwasawa et al. 



2Ullb[ ). All three nuclei have optical peaks in with the 



W nucleus being 2x as bright as th e E nucleus in 2 kpc 
square apertures ( Armus et al.|199 Q'), and the sum of the 
W+SE nuclei being 4 x as bright overall as the E galaxy 
( Murphy et al.|p"996 ). The W nucleus also contains 95% 
of the system's tota l Ha flux (uncorrected for extinction; 



Armus et al, 
are together 



1990|). In the K band, the W+SE nuclei 
6x as bright as the E nucleus (based on 



2MASS K20 magnitudes retrieved via NED), and the W 
nucleus i s in turn sign ificantly brighter than the SE nu- 
cleus (Murphy et al.||l996 ). 

Judging from Infrared Space Observatory images, the 
E nucleus is very faint compared t o the W+SE complex 
at 7 and 15/im ( [Dale et al.]|2QQQ ). Most of this MIR 
flux in turn arises from the W nucleus, as seen in Spitzer 
archival images at 4— 24 /im. The W nucleus is the source 
of the prodigious FIR emission from this system, which 
is obvious from archival Spitzer 70 and 160/im images. 

The E nucleus is a red spheroidal system with no star 
formation, as evident from HST optical images and the 
archival SDSS spectrum, which show a smooth light dis- 
tribution and an absorption-line spectrum. The exact 
nature of the SE nucleus is unknown, but its blue color 
is apparent in Figure [TT] The W nucleus is classified op- 
ticall y as H Il/composite (jVeilleux et al.|1995 Yuan et al. 
2Q1Q[), but it shows no hign-excitation IVlIK emissio n lines 



I^arr ah et al.^2QQ 7|) or X-ra y signatures of an AGN ( Teng 



■Veilleux 2010; Iwasawa et al. 2011bJ. Based on FAH 



strength and MIR spectral shape, , |Vehleux et al.| ([2009) 
determine an AGN fraction of 0.17! 

We conclude that a starburst in the W nucleus of 
F10565+2448 is producing most of the energy in the sys- 
tem and that an AGN, if any is present, contributes very 
little to this energy output. 

4.2.2. Kinematics 



As discussed in Shih & Rupke (2010), the ionized gas 
rotates in the same sense as the molecular gas, in a dis k 
of PA - 100° and i - 20° ( |Downes fc Solomon] [19981 ). 
The near side of the disk is not obvious because of the 
irregular morphology of the galaxy, but, as we discuss 
below, it is likely to be the N side. 
A neutral outflow was discovered in this system via nu- 



y 

clear spectra ( Heckman et al.|[2000 Martin|[2005 Rupke 
et al. 2QQ5c[)^ and i t is visib le m molecular gas, as well 
(|Sturm et al. 2011). Martin (2006) showed that this neu- 
tral outflow extends eastw ar d at l east 8 kpc from the 
nucleus, while jShih fc Rupke (2010) showed that its pro- 
jected extent is at least 4 kpc from the nucleus i n most 



other directions, as well. Shih & Rupke (2010) found et al. 



that the neutral gas outflow is dusty and has a velocity 



They concluded that the outflow may be decelerated in 
the E due to ambient interstellar gas, since significantly 
higher velocities are observed W of the nucleus. 

We verify the extent and strongly asymmetric velocity 
structure of the neutral and ionized gas found by |Shih 



fc Rupke (2010). The neutral gas velocity varies slowly 
with projected galactocentric radius, as we show in Fig 



ure 



16 



At a given position angle (E of N with respect 
to tKe" galaxy center), these figures show either flat ve- 
locity profiles or perhaps a modest decline in ^50% with 
projected radius. SRFW models are able to reproduce 
this behavior as long as the velocity is quite different 
E and W of the nucleus and the wind radius is large 
enough (at least 5 kpc). The radial profiles of Vgs% do 
show a noticeable decline with radius. If both ^50% and 
Vgs% represent bulk flow of gas at similar locations in the 
wind, then this suggests that the gas radius is also not 
too large (^10 kpc). Alternatively, the changes in Vgs% 
may reflect the decline in unresolved bulk or turbulent 
motions with radius. 



Compared to the analysis of Shih fc Rupke (2010), the 
current study increases the spatial information available 
for this source by Nyquist sampling the seeing disk. With 
0'!3 spaxels, a minor axis outflow on scales of ~1 kpc is 
evident in the data (Figure 17). The primary evidence 



for this is in the morphology and extent of the largest 
FWHM, highest velocity regions in the neutral and ion- 
ized gas. If we assume that the blue clumps to the N 
of the nucleus represent unobscured star forming regions 
in the near side of a galaxy disk, then the high velocity 
regions extend directly S along the minor axis of the CO 
nuclear disk. We would expect blueshifted emission on 
the S side of the disk, as is observed. The angle that 
the wind makes with respect to the line of s ight should 



be nea r that of the disk's inclination (~20°; Downes fc 



Solomon 1998). The minor-axis flow must then be rela- 
tively collimated (total opening angle, from one edge of 
the wind to the other, of < 2 x 20°), since we observe it 
primarily in projection away from the line of sight. 

The regions of highest velocity are strikingly similar 
in morphology to the dustiest regions of the disk, if we 
assume that the HST color just S of the nucleus is due 
to variable obscuration (consistent with the patchy and 
filamentary structures of the colors in this region). The 
largest values of FWHM in the gas (up to 800 km s~^) are 
coincident with the dustiest region of the galaxy (and re- 
gion of highest Na I D equivalent width). Dust filaments 
extend S along the minor axis, coincident with regions of 
high velocity, suggesting that dust is being entrained in 
the wind along this direction. 

4.3. Mrk 231 



Mrk 231 is a late stage merger in which the nu- 
clei have al ready coalesced ( |Surace et aL]|1998[ |Veilleux | 
et al. 2002[ ). It is the nearest CjSO ( [Bofesenberg et al. 
1977p , as well as the nearest broa d absorption line (BAL) 



( [Adams fc Weedman|[1972| . The kpc-sc ale neutral 
outflow in this system was discovered by Rupke e t al.| 
(2005a), a nd was more rec ently found to extend in all 
directions ( Rupke fc Veilleu x 2011) and have a signifi- 
cant molec ular component ([Fischer et al.|2010a| Feruglio 



Mr 



2010[|Aalto et al.[[2QT2{ [Cicone et aT 



i 231 is dominated by a central point source (3 x 



profile inconsistent with a symmetric minor-axis outflow. brighter than the host at H] ^Veilleux et al.,^2006J that 
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has a QSQ spectrum wi th broad B aimer and Fe lines 
( Boksenber^ et al. 1977|). It al so hosts a massive CO 
disk ( [Downes &: Solomon 1998) with a high star for- 
mat iorrTaIe7noaieT"oirTQe presence of PAH emission, 
[Ne II] 12.81 jam emis sion beyond that r equired f or an 



AGN , and MIR colors ( [Genzel et al.|1998 
[2QQ9| ). The AGN produces 



Veilleux et al. 

7U% of the galaxy's lumi- 

nosity ([Veilleux et al.||2QQ9|) . 

The kinematics of the disk rotation and outflow in 



Mrk 231 are discussed in'Rupke & Ve illeu^ ( [2Qll| ). We 
present the entirety of the Mrk 231 data set nere tor easy 
comparison with the rest of our sample, but do not dis- 
cuss it further. Figures 18-21 show maps and velocity 
distributions for this system. 

4.4. Mrk 273 

Mrk 273 is a late merger with a binary nucleus and an 
AGN. The bolometric luminosity of the galaxy is proba- 
bly dominated by a starburst, but the AGN contribution 
is still signiflcant. The kinematics in this system are com- 
plicated, but clearly show a high-velocity outflow. Maps 
of derived properties and velocity distributions can be 
found in Figures [22]-[26) 



4.4.1. Power Source 

The nuclear properties of Mrk 273 are enigmatic. Two 
nuclei (typically labeled N and SW ) are e vident in the 
NIR and M IR (Majewski et al. 1993; Knapen et al. 1997j 
Scoville et al. 2000; Soifer et al. 2000; see Figure 1231) 



I'he N nucleus contains most oi the molecular gas m 
the syst em within an extreniely co mpact core (radius < 



I II t ^ 

120 pc; [Downes fc Solomon[[1 998 L and is the s ource of 

most of the MIK luminosity ( [!Soiter et al.[[2000[ ). In the 

radio, Mrk 273 resolves into a dittuse radio source with 



two primary peaks ( Carilli fc Taylor|[2000 ) . The N radio 
peak (coincident with the N nucleus) is further resolved 
as multiple sources at the highest resolutions, and the SE 
radio peak is i nconsistent with an AGN based on spectral 
slope dBondi et al. 2005). The optical line emissio n from 



the center of Mrk 273 classifies it as a Se yfert 2 ([Koski 
[1978 j Veilleux et al.[jl999aj . Using IFS, jUolina et al ' 



locate the Seyfert 2 emission near, but to the SW 

ot, the two nuclei. 

Recently, the SW nucleus has been identified as an 
AGN through localiza tion of the hard X-ray point source 



( [Iwasawa et al.|2 011b| ) , which is C ompton thick bas ed on 
spectral modeling trom 2—40 keV (Teng et al. 2009). The 



AGN shows strong hi^h-ionization narro w-line emission 



1999b) and stroni 



Armus et al 



n|2007l 



in the NIR and MIR: possible [Si VI] ( Veilleux et al. 
' ^ [Ne V] and [O IV] lU 

ial[2on 



Genzel et al. 
Ml. The ra- 



Veilleux et 

tios of these high-ionization lines to a low-ionization line, 
[Ne II], point to a do minant contribution of the AGN 
to the FIR luminosity (Veilleux et al. 2009). Combined 



PAH and emission-line diagnostics su ggest that the star- 
burst and AGN contribute similarly ( Genzel et al.[ [1998| 
[Veilleux et al. 2009), while MI R continuum d iagnostics 
point to a dominant starburst ( Veilleux et al.[[2009 ). 

NIR IFS of the inner kiloparsec ol Mrk 273 local- 
izes high-ionization coronal line emissio n ([Si VI]) to th e 



SW nucleus and the SE radio source (U et al. 



2012|). 

larbor 



These new data suggest that the N nu cleus may ^ 

a heavily-obscured AGN (see also Iwasawa et ar][2QTTb ) 



that photoionizes the gas in the SE radio source. Alter- 
natively, high- velocity shocks may produce this emission. 
Regardless, the existence of an obscured AGN in the N 
nucleus remains to be investigated. 

The concentration of CO and MIR continuum emission 
in the northern nucleus would seem to imply that the 
SW AGN contributes little to the bolometric luminosity 
of this dusty system. However, MIR dianostic s suggest 



otherwise, with an average AGN fraction of 0.46 (Veilleux 
|et al.|2009| ). This conclusion relies on the combined MIK 
spectrum of both nuclei, and rests primarily on the high 
emission-line ratios. The bolometric AGN luminosity 
inferred this way is entirely consistent with the X-ray 
emission of the SW nucleus if we compare it to optically- 
selected QSOs. If the AGN luminosity were lower by 
more than a f actor of 2 — 3, then L2 -10 keV /Lagn would 
be too high QTeng fc Veilleux[[2QTQ ) . 



4.4.2. Basic Kinematics 

The gas kinematics in Mrk 273 are the most complex 
in our sample. The atomic and molecular gas kinemat- 
ics centered around the N nucleus are characterized at 
100 pc to kpc scales by a rotating disk of gas ( Do wnes] 
fc Solom on|1998[ [Cole et al . 1999; Carilh fc Tayl or 20W 
Yates et al. 20^. The position angle of the line of nodes 
ol this disk is 7(7—90°, an d it is inclined by 45 — 55° with 



respect to the line of sight (Downes fc Solomon 1998 Cole 
et al.[[1999| [U et al.[[2012[ ). I'he disk is rotating rapidly 



with a projected velocity gradient of 400 — 500 km s~^ 
over a few hundred pc (velocity increasing from east to 
west; [Downes fc SolomonI [19981 [CoITet al.| |1999[ [Cax^ 



iUi fc Taylor[[2000| [Yates et al [[2000p . On larger scales" 
there is a molecular feature that extends 5 kpc N and 
one that extends 5 kpc S of the N nucleus; the N feature 
is redshifted by ^3 00 km s~^. and the S feature is near 
systemic (Downes fc Solomon[p!998J . 

There are signiticant departures from simple E-W ro- 
tation, from sub-kpc to kpc scales. The spatially re- 
solved CO lines a re broad, and the in tegrated profile 
spans 1000 km s~^ ([Downes fc Solomon 1998| ). T wo com 



ponents eme rge in the [O III] 5007 A spectra of Colina 



et al, 



(1999), and they decompose these into a systemic 



component and a component that is redshifted in the E 
and blueshifted in the W, with a peak-to-peak amplitude 
of 2400 km s~^. They suggest that the latter component 
represents a starburst-driven superwind. 

The present IFS data represent a significant improve- 
ment in sensit ivity, spatial resolu tion, and spectral reso- 



lution over the Colina et al. ( 1999| dataset, as well adding 
information on the neutral gas phase. 

It should be noted that the CO/H I gas disk is tiny 
compared to the F OV of our map s, subtendin g only 1 — 



2 CMOS sp axels ( jDownes fc Solomon,, 1998 Carilli 
Taylor[[2000[). 



i'he ionized gas profiles are extremely broad in many 
spaxels, and are in most cases best fit by two or three 
components. Fig. [27] shows several examples of these 
fits, and Fig. [28| shows the number of components in each 
spaxel. We omer these components into one of four cat- 
egories as follows: (1) rotating, (2) broad and redshifted, 
(3) broad and blueshifted, and (4) isolated clouds with 
narrow profiles. In a two-component spaxel, we assume 
the narrower component is rotating, though we tweak 
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some spaxels by hand because the FWHM of the two 
components are similar. The second component is then 
assigned to the broad and redshifted bin or broad and 
blueshifted bin depending on its velocity with respect to 
the rotating component. The isolated clouds with nar- 
row line profiles are easily distinguished from the rotat- 
ing narrow component because of the significant red- or 
blueshifts of the isolated clouds, and from broad com- 
ponents by linewidth and spatial footprint. If there are 
three components in a spaxel, we generally choose the 
central component to be rotating, though there are some 
exceptions. 

The rotating i onized component follows the CO (1-0) 
rotation curve of Downes & Solomon (1998) very closely, 
including the twisting of the systemic isovelocity contour, 
the peak redshifted velocity just W of the nuclei, and 
the redshifted gas to the N. The rotating component is 
on average narrow compared to the other components 
((FWHM) = 300 km s-^). The line ratios are low along 
a NNE/SSW axis, with log([N II] /Ha) - -0.6 - and 
log([0 I] /Ha) in the range -1.8 to -0.8. These ratios 
increase to 0.3 and -0.5 to the ESE and WNW. 

The redshifted and blueshifted b road, ionized compo - 



ttii 



1999): 



nents are broadly consistent with [Colina et al. 
the blueshifted components lie primarily W of 
clei, while the redshifted components lie to the E. How- 
ever, our higher spatial and spectral resolution shows 
significant substructure. Projected velocities rise away 
from the nucleus, peaking at a remarkable v^q% = 

— 1500 km s~^, 4 — 5 kpc NW of the nuclei. However, the 
broadest components (FWHM up to 1500 km s~^) with 
the highest [N II] /Ha and [O I] /Ha line ratios (indicative 
of shock ionization) arise directly N and S of the nuclei. 
(These locations also display modest peaks in \v^q%\ of 

- 1000 km s-^) 

The neutral gas in Mrk 273, as probed by Na I D, is 
backlit by the stellar continuum across most of the N-S 
extent of the CMOS FOV in a band that is 1-2 kpc wide. 
There are few departures from velocities near systemic. 
However, there is blueshifted gas north of the nuclei. In 
the SRFW model calculations for this system, we include 
only spaxels in which the neutral gas is blueshifted from 
the ionized gas narrow component in that spaxel by more 
than 50 km s~^. 

The blueshifted ionized and neutral gas in the N/NW 
of the system is clearly outfiowing, as is the red- 
shifted ionized gas in the S/SE. The velocities (reaching 
1500 km s~^ in projection) are too high to be explained 
simply by tidal motions, and are much larger than the 
projected velocities in the narrow component (between 
—350 and 350 km s~^). However, the redshifted and 
blueshifted ionized and neutral components in other lo- 
cations have more modest velocities ('^50% and FWHM 
of a few hundred km s~^) that could in principle be ex- 
plained by tidal motions. 

The eviden ce for a bip o lar ou tfiow aligned E-W, as 
proposed by ^Colina et al. (1999), is not the most likely 
explanation based on the current data, though a varia- 
tion on this theme emerges from our analysis. The data 
conclusively point to a bipolar nuclear superbubble in the 
center of the system, aligned N-S. Gas further from the 
nucleus, if it is indeed outfiowing, then traces gas blown 
out from this bubble or earlier ejections from the nuclei. 



4.4.3. Nuclear Superbubble 

Position-velocity (PV) diagrams of the line emission 
in the center of Mrk 273 reveal the classic signature of 
a bipola r supe rbubble (Figure 29). Applying our BSB 
model (§3.4.2) and optimizing its parameters yields the 
values given m Table [5] The model is overplotted in the 
PV diagrams, and the two-dimensional projected velocity 
fields and bubble shapes are compared to the data in 
Figure 30 The bubble is fairly symmetric in emission 

The absorption lines in Mrk 273 are not of high enough 
signal-to-noise near the possible superbubble to accu- 
rately parameterize its neutral gas content. There are 
higher absorption line velocities near the edges of the 
near-side bubble in the N, suggesting that some neutral 
material is entrained by an expanding bubble. 

Overall, a bubble shape fits the inner radii of the wind 
well (at projected radii < 2 kpc). However, the high 
velocities seen immediately adjacent to it but at larger 
radii in the NW and SE strongly suggest that the bubble 
is broken on its surface and gas is escaping the bubble in 
some directions. This gas has as high or higher veloci- 
ties than in the bubble itself, perhaps indicating further 
acceleration by a hot medium that could be blowing the 
bubble. 

There are also areas of high velocity directly E and W 
of the nucleus, along a position angle perpendicular to 
the bubble axis. These may represent earlier episodes of 
gas ejection, or ejection due to an unrelated process. 

The fitted outfiow center is formally located between 
the two nuclei, 0'.'4 (~300 pc) from the N nucleus and 0'!7 
(~500 pc) from the SW nucleus. However, the errors in 
the bubble center are 1 spaxel (0^3), so an origin at one 
or the other of the two nuclei is plausible. If in fact the 
outfiow arises between the two nuclei, then it could be 
powered by a combination of energy from the AGN and 
starburst. However, the outfiow may also be asymmetric 
(extending further on the near or far side due to, e.g., 
asymmetric acceleration), lending further uncertainty to 
the true outfiow center. 

If the bubbles were centered on the N nucleus, their po- 
sition angle (—5°) would be consistent with collimation 
by the tilted CO/H I disk at that location ( PA 70 - 90°; 



Downes fc Solomon][T998[ |Cole et al.||l999| ). The disk's 
inclination is i = 45" from CO modeling (Downe s fc| 
SolomonI 19981), which is inconsistent with the superbub- 
ble inclination of 75±5°. However, a higher disk inclina- 
tion is not out of the q uestion, given the large velocity 
gradients in the nucleus (Cole et al.||1999 Carilli & Tay- 



lor 



2000 Yates et al. |2000p . A more recent analysis of 
NIkTFS data, taken with adaptive optics, finds a more 
inchned disk {i ^ 55°) based on ionized gas lines (U et al. 
2012). 

The fitted outfiow center is further from the SW nu- 
cleus. If the bubbles were centered on the SW nucleus, 
it is unclear what would collimate the outfiow. Further- 
more, the N disk has a large gas reservoir to couple with 
the energy driving the outfiow, while the SE nucleus does 



not ([Do wne s fc Solomon|1998[ |Cole et al.|1999[ [Carilh 



T^or||2000| ). 



Our structural analysis is inconclusive on the power 
source of the entire outfiow (superbubbles + extended 
fiow), but more consistent with the superbubbles being 



10 



RuPKE & Veilleux 



powered by the N nucleus than the SW nucleus. Data 
of higher spatial resolution are needed, as well as con- 
straints on the presence of an obscured AGN in the N 
nucleus. As we discuss below, how ever, the high veloci- 
ties suggest an AGN-driven flow (§5.4.3). 



4.5. VV 705 

VV 705 is a binary merger with nuclear separation of 
7.5 kpc, based on inspection of an archival HST Wide 
Field Camera 3 F160W image. Maps of derived proper- 
ties and velocity distributions for each nucleus are shown 



in Figures 31-40 The superbubble emerging from the 



NW nucleus is one of two in our sample. 

4.5.1. Power Source 

The NW nucleus contains most of the luminosity in 
this system at multiple wavelengths: 8.44 GHz (NW:SE 
luminosit y ratio of 3:1; [Condon et al.| |1991), 2 — 10 keV 
(8:1; Iwasawa et al.||2011a|), and 4 — 24 /im (archival 
SpitzerVKA^ and MIPS 24/im data). Though exact mea- 
surements are not available, at least 90% (and perhaps 
closer to '^99%) of the 24 fim emission emerges from 
the NW nucleus, based on the strength of the 24 /im 
PSF of the NW nucleus. This strongly suggests that 
this nucleus is the source of most of the FIR luminos- 
ity of the galaxy, as well. This would be consistent 
with it hosting most of the molecula r gas (~10^^ Mq; 
Sanders et al.|199T Chini et al.ll992|), based on the nar- 



row FWHM of CO (100 km s 
the measured CO velocity n ear t 



[Chini 

hat of 



et al. 



theJN 



1992) and 



W nucleus 



{cz = 12070 - 12090 km s"^ ; Chini et al.,1992 Yao et al. 



2003[ and this work). 

Both nuclei are classified as compo site by [Yuan et al.] 
(12010 ), based on spectra obtained iDy Kim et al. ( |1995l T^ 
'Fhis IS confirmed for the NW nuc leus by nuclear data 

( 2006 J . However, there is 
ler nucleus; X-ray i maging 



from Moustakas & Kennicutt 
no evidence for an AGN in eit 

shows soft spectra in both nuclei ( Iwasawa et al . |[2Q 1 la), 
MIR spectra reveal no [Ne V] or [O IVJ tine structure 
lines ( Dudik et al.||2009l ), and comp act radio emission is 
consistent with a nuclear starburst (Smith et al. 1998). 

In the absence of evidence for an AGN, we conclude 
that any outfiows in this system will be powered by 
star formation. We assume an AGN fraction of 10% for 
consistency with the other two galaxies in our sample 
without luminous AGN, but the actual fraction could be 
lower. 

4.5.2. Kinematics 

The SE nucleus shows ordered rotation with non- 
circular motions superimposed. We estimate a line of 
nodes PA of —30°. The peak-to-peak velocity amplitude 
along this PA reaches 250—300 km s~^. The near side 
of the disk is likely to the W, given that this side of the 
nucleus is less obscured (Figure 36). The line of nodes 
and rotation amplitude of the NW nucleus are more dif- 
ficult to quantify due to tidal motions and a near face-on 
orientation. We do not attempt to estimate a PA or in- 
clination for this nucleus. 

Non-circular motions and multiple narrow components 
in the SE nucleus include a region of LOS overlap be- 
tween the two gas disks, 3 kpc N of the SE nucleus. 
Multiple components are also evident within 1 kpc of 



the SE nucleus; some of these components have FWHM 
exceeding 600 km s"-*^. These broad components have 
velocities within 100 — 200 km s~^ of the narrow compo- 
nent in the same spaxel. They may arise in a weak out- 
fiow along the minor axis and/or with infalling shocked 
tidal debris. Either of these interpretations are consis- 
tent with the observed elevated FWHM and line ratios, 
but the data is not sufficient to distinguish between these 
possibilities. Regardless, any outfiow in this nucleus is 
not powered by the primary source of luminosity in the 
merger. 

The NW nucleus shows evidence for a superbubble 
emerging almost directly along the line of sight. This 
bubble is most evident in absorption, as seen in t he PV 
diagram in Figure 41 Applying our BSB model (§|3.4.2) 



and optimizing its parameters yields the values given m 
Table [H The fits and residuals in two dimensions are 
shown in Figure[42] 

The superbubble model fits the inner 2 kpc well, but 
outfiowing gas at lower velocities is also visible NE of 
the bubble and, to a lesser extent, to its SW. The best- 
fit parameters show that the outfiow emerges almost di- 
rectly along the line of sight and reaches a velocity of 
'^50% = —500 to —600 km s~^ along the bubble axis. 
The elongation of the bubble is fairly uncertain, and is 
consistent with either a symmetric bubble (~2 kpc in di- 
ameter) or one that is longer perpendicular to the galaxy 
disk by factors of 2 — 3. 



The outfiow is also seen in emission (Figure 33). How- 
ever, the ionized gas only shows acceleration m a frac- 
tion of the outfiow area that is seen in absorption. The 
broad component in most spaxels has the same veloc- 
ity as the narrow, rotating component. As a result, in 
the SRFW model calculations for this system, we include 
only spaxels in which the broad component is blueshifted 
from the narrow component in that spaxel by more than 
50 km s~^. 

4.6. F17207-0014 

This system is a coalesced merger with a diffuse, patchy 
optical morphology. We present maps and velocity dis- 
tributions based on our GMOS data in Figures [43] 



47 



Fl 7207— 0014 is the only system in our sample without 
a detected ionized outffow, though it does reveal a large 
scale neutral gas outffow. 

4.6.1. Power Source 

All of the data on this system is consistent with an 
obscured starburst. The nu cleus is obscured in the opti 



cal but reso lved in the NIR (| Murphy et al.|1996 



'et al.^"2000; Figure [44l). It is a diffuse source in t 



Scoville 
le radio 



and MIR (Soifer et^aT 2000; Momjian et al.|2003l Baan & 
Klockner 2006) , no |Ne V| or |0 IVJ emission is detecte( 
in the MIR (Farrah et al 1200 71), the AGN fraction based 



3 



on MIR diagnostics is only 0.11 ( [Veilleux et aI||2009D 
and the X- ray data do not show an obscured AGN (e.g. 



Francesch ini et al.|[2QQ3l |Teng &: VeilIeiS[[2QTQl [Iwasawa 
|et al.||2011b^ 



4.6.2. Kinematics 

The CO and H I disk in F17207-0014 subtends 1 
along a line of nodes position angle o f ^120° (Downes 
'fc Solomon| [1998[ [Momjian et al.[ [2003[ ). [Arribas"SrCQF 



ina^ (2003) presented ionized gas maps of this system and 
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showed that the molecular and neutral kinematics are 
seen in the ionized gas, as well. The rotation of the nu- 
clear ionized disk is roughly centered on the NIR peak, 
but at larger scal es the isovelo city conto urs twist away 
from this pattern ( Arribas fc Colina|2QQ3 ). The kinemat- 
ics on ^al actocentric radii ^,1 kp c probably reflect tidal 
motions ( [Arribas fc Colina 2003). 

A neutral outflow in this system was discovered by 
flupke et all ( |2005c| ) and Martin (20 Q5|, and shown t o 
extend over 5-6 kpc ( |Rupke et al |20d5c| partin|20 06') . 



Martin ( 2006 ) found no variation in the outffow velocity 
(with respect to the rotating ion ized ^as) along a posi- 
tion angle of 170°. Recently, Westmoquette et al. (2012) 
presented ionized gas maps of F17207— UU14 that showed 
two velocity components in the ionized gas in the center 
of this system. However, the two components have very 
similar velocity, and are not obviously an ionized outflow. 

Our new maps of the ion ized and neutral gas in 
F17 207-0014 improve on the [Arribas fc Colina|^2003 ) 
and Westmoquette et al. ( 2012p studies by adding the 
neutral gas kinematics and increasing the spatial resolu- 
tion. Our ionized gas maps indeed resemble higher reso- 
lution versions of the Arribas fc Colina (2003) maps. We 
do not detect two ioni zed gas components, as in West- 
moquette et al. ( 2012|), but th e much larger spaxels in 
the I Westmoquette et al. (2012) data may increase their 
sensitivity. 

We detect extended neutral gas that is everywhere out- 
flowing in this system, consistent with previous nuclear 
and extended detections. However, we show that this 
outflow extends in all directions, and to galactocentric 
radii of at least 4 kpc (with a maximum extent of 6 kpc, 
as in Martin 2006|^ 

Contrary to Martin (2006), we find that the outflow 
velocity increases with increasing projected radius (Fig- 
ure 48). This is inconsistent with our fiducial SRFW 



model, which predicts decreasing velocity with increas- 
ing projected radius. An alternative scenario is that the 
projected galactocentric radii represent the true galac- 
tocentric radii, and that the outflow speed is increasing 
with radius. 

5. SAMPLE TRENDS 

Our sample of six merging systems reveals an extended 
outflow in every case. In flve systems there is an ionized 
outflow (F17207-0014 is the odd galaxy out), and four 
of these are new ionized outflow detections (th e Mrk 273 
outflow was previously known; [Colina et al.| fl999). In 
six systems there is a neutral outflow. 'I'hree oi these are 
new detections (F08572+3915, Mrk 273, and VV 705), 
despite two of these three having been observed in neutral 
outflow surveys using single-aperture spectroscopy. 

In this section we combine the results from our six- 
galaxy sample to draw broader conclusions about the 
multiphase structure and power sources of galactic winds 
in ULIRGs. These conclusions are summarized in §[6] and 
rely in part on the analysis in Figures 49 - 56 



5.1. Wind Structure 
5.1.1. Ionized Gas 

We know a great deal about the structure of winds in 
local disk galaxies hosting starburst-driven superwinds 
( Veilleux et al.^ 2005} . These winds typically emerge 



along the minor axis of the disk as either a bipolar su- 
perbubble or a bicone. The ionized gas in such BSBs 
or bicones is often limb-brightened and fllamentary, sug- 
gesting it lies on the edges of the bubble or cone (e.g. 



■Heckman et al.n9871 [Bland fc T ully 19881 [WTleux et al 

1994; Martin 1998, among many examples) . In the model 
of a bipolar wind, these bubbles or cones have opening 
angles that increase with radi us, reaching angles (o n one 
side of the disk) of 45 - 100° ( [Veilleux et al.[ [2005). 



The structure of ionized gas outflows in local mergers 
is less understood, though the two nearest and most- 
studied examples have possible minor-ax is and/or bipo- 
lar flows: Arp 220 ([Heckman et al.|1987t and NGC 6240 
dVeilleux et aL]|2003r[Bush et al.||2005| r The kinematics 
of the ionized gas m mergers are made complex by the 
preval ence of multiple overlapping disks and tidal mo- 
tions (Colina et al. 2005), as well as outflows. However, 
careful spectral decomposition of the emission li nes is 
able to separate gas rotation from other motions ([Shih 
Rupke[|2QTQl [Rupke fc Veilleux|[2QTT| [Westmoquette 



Kupke||:^UlUt |Kupke ^ Veilleux||:^Ull[ [Westmoquette 
et al. f2Q12[ and this study). The narrow component in 



the ionized gas (always the strongest in flux overall) fol- 
lows the rotation of the resolved kpc-scale CO disk (when 
CO has been previously observed via inter ferometry at 
arcsecond resolution, as in 4/6 systems in our sample). 
The good match between the ionized and molecular gas 
in the inner 1 — 2 kpc of these systems is somewhat sur- 
prising, sin ce these disks ar e dusty and probably opti- 
cally thick (Downes fc Solomon 1998); thus, the ionized 
gas is likely to trace the outer "skin" of the disk rather 
than its core, or emerges from a patchy and inhomoge- 
neous medium. Where the molecular gas is more ex- 
tended (e.g., Mrk 273), there continues to be a good cor- 
respondence between the bulk motions of the molecular 
and ionized gas. 

The linewidths in the narrow component of the ion- 
ized gas are relatively large (the average FWHM in this 
component in each of the eight nuclei in our sample 
is in the range 100 — 300 km s~^, and typically near 
200 km s~^; see, e.g.. Figure j5|. Similar linewidths 
are seen in the resolved m olecular gas in these systems 
( Downes fc Solomon||1998 ). These large linewidths may 
result from several effects. They do not likely result 
from beam smearing, since the narrow component ve- 
locity flelds typically have smaller amplitudes than this. 
Sometimes overlapping disks and/or regions where two 
disks are colliding (as in F08572+3915 and VV 705:NW) 
amplify the linewidths in isolated regions as a result of 
the combined motion of the two disks or the tidal motions 
of collision. In other systems, the linewidth may be in- 
flated by either energy injected by star formation and/or 
AGN activity in the inner disks, or on larger scales by 
tidal motions resulting from the ongoing merger. 

On scales larger than the central 1 — 2 kpc, the ionized 
gas continues the rotational pattern of the inner disks, 
but often with isovelocity twists superimposed or changes 
in velocity amplitude. It is likely that the motions of the 
ionized gas on large scales in most of these systems bet- 
ter reflect the ongoing merger of the two initial systems 
rather than the inner disk of either nuclei or the nuclear 
disk of the coalesced system. The ionized gas motions 
at projected galactocentric radii of 3 — 5 kpc thus result 
from a combination of the internal angular momenta of 
the original two systems with that of the merger orbit. 
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The interpretation of the broader and blueshifted or 
redshifted velocity components in the ionized gas is rel- 
atively straightforward. In emission and without other 
information, outflow and inflow are in principle indistin- 
guishable. However, as we have shown for Mrk 273 and 
VV 705, outflows can be distinguished from inflow via 
structural modeling, since BSBs are obviously present. 
A multiphase, multi-kpc minor axis outflow is also likely 
present in F 10565+2448, based on the morphological 
connections between the ionized gas kinematics, neutral 
gas kinematics, and nuclear dust filaments. In the cases 
of F08572+3915, the large-scale kinematics in Mrk 273, 
and the inner kpc of Mrk 231, the observed velocities 
on these scales are simply too high to be explained by 
t i dally- induced inf low . 

|Soto et al. (2012) also argue that, in ULIRGs, the ob- 
served ionized gas must be on the near side of the galaxy, 
and thus outflowing if blueshifted, since far-side ionized 
gas would be obscured by dust. There is some evidence 
from our data against this conclusion, but there is also 
evidence to support it; the reality may be more nuanced. 
First, the continuum obscuration by dust (as revealed 
by multicolor iJ^T imaging; e.g.. Figure Il7|) is patchy. 
Nonetheless, optically thick molecular disKs (not traced 
by the optical continuum) reside at the nuclei of these 
systems and probably obscure any far side ionized gas 
very near the nucleus. Second, we do observe the reced- 
ing side of the superbubble in Mrk 273. However, this 
bubble is highly inclined and extended; in most other 
cases where we see evidence consistent with a near-side 
minor axis ionized flow, evidence for the far-side flow is 
weaker or absent (see discussion below and of individual 
objects above). The good corres pon dence between ion- 
ized and neutral gas velocities (§5.3) also suggests that 
the ionized gas is primariliy located on the near side of 
these systems, where the absorbing gas must lie. Over- 
all, it is clear that spatially resolved observations of these 
systems, as we present here, are crucial to resolve these 
structural degeneracies. 

The ionized outflows we observe are in all cases consis- 
tent with being collimated by a molecular disk within the 
galaxy, emerging along the minor axis of that disk. In 
the three cases where CO disks and ionized outflows are 
observed (F10565+2448, Mrk 231, and Mrk 273), there 
are outflows reaching 1 — 2 kpc away from the nucleus 
that have directions and inclinations consistent with the 
CO minor axis. In the other two cases of ionized outflows 
where CO disks are present (but have not been observed 
via interferometry ) , the direction and inclination of the 
outflows are consistent with the disk minor axis as in- 
ferred from the rotation of the ionized gas. 

In particular, in the two cases where we fit the 3D 
shape of the outflow, the kinematics are consistent with 
2 kpc diameter BSBs, which would formally have (one- 
sided) opening angles of 180°. In reality, the footprint of 
the outflow is not a point, such that the opening angle 
is smaller than this. For Mrk 273, let us assume the 
CO core/disk edges from Downes & Solomon' fl998) as 
the extent of the launching region (radii 100 — 400 kpc, 
depending on the measure) . If we draw a line connecting 
the CO disk edge with the bubble edge at a disk height of 
1 kpc, then we infer an opening angle of 30 — 40°, similar 
to local superbubbles. 

Outflowing ionized gas is also observed away from the 



minor axis, particularly in the cases of Mrk 273 and 
F 10565 +2448. The morphology of these extended flows 
is more ambiguous, though in Mrk 2 73 it may in part 
be breakout from the superbubble (§4.4.3). It is pos- 
sible that these are traces of earlier bubbles or ejection 
episodes that have moved to larger radii. 

This idea is buttressed by the timescale for wind evo- 
lution over the galactocentric radii probed by our ob- 
servations, which is tdyn ^ Rwind/v50% = 10 Myr (for 
Rwind = 3 kpc and '^50% = 300 km s~^). The timescale 
for the orbital evolutio n of the galaxies in our sam ple is of 
order ^,100 Myr (e .g., Barnes & Hernquist 1996( [Mihos 
Hernquist 1996); i.e., an order of magnitude larger. 
I'hus, if we assume that we are observing the six galax- 
ies in our sample at a random time late in the merger, 
then the existence of minor axis flows at roughly similar 
scales (1 — 2 kpc projected radii) in each system suggests 
that they are a relatively continuous phenomenon, either 
in the form of multiple ejection episodes or a continuous 
flow. This in turn further suggests that the removal of 
gas from the nuclei of mergers is an efficient process. We 



return to this subject in § 5.5 



5.1.2. Neutral Gas 

The prevalence of minor- axis outflows collimated by 
a CO disk is seen in neutral, as well as ionized, gas. 
F10565+2448, Mrk 231, Mrk 273, and VV 705:NW ah 
show evidence of expanding neutral flows along the minor 
axis, roughly consistent with the morphology of the ion- 
ized flows. In the neutral gas, however, there are excep- 
tions: F08572+3915 and F17207-0014 do not obviously 
show flows along the CO minor axis, despite revealing 
outflowing neutral gas at other locations. Furthermore, 
in F10565+2448 and VV 705:NW, there are extensive 
areas of outflowing neutral gas away from the minor axis 
flow. 

In the case of neutral gas, there is structural ambi- 
guity for two reasons: (1) the neutral gas requires a 
strong continuum background source, so the absorption- 
line signal-to-noise ratio decreases away from the galaxy 
nucleus; and (2) the Na I D resonant line may suffer 
from emission-line filling near systemic ( Prochaska et al. 
2011 ). We try to minimize the former by deep exposures, 
but there is still a limiting radius at which the signal 
drops below detectability. Though we cannot provide 
quantitative constraints, we believe the latter to be min- 
imal, since the strong dust obscuration in these systems 
would prevent outflowing emission from the fa r side of 
the outflow from reaching us (see discussion in 



5.1.1 for 



more on far-side line emission). This qualitative conclu- 
sion is supported by the lack of any observed redshifted 
Na I D emission in these systems, even at large radii 
where the dust obscuration is smaller than in the nuclei. 

5.1.3. Maximum Extent and Opening Angle 

The extended outflows in our dataset are consistently 
large, stretching across much of the FOV. Five of the 
six systems in our sample show outflowing gas up to the 
edge of the FOV, at scales ranging from 3 — 5 kpc from 
the nucleus. Thus, the sizes inferred from our data are 
ultimately lower limits to the overall wind sizes in these 
systems. As shown by "Martin" (2006") , outflows continue 
to larger radii in long slit spectra. Using larger radii 
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in our SRFW model would only increase the measured 
masses, momenta, and energies. 

The location of the neutral gas in nearby galaxies is not 
fully constrained because it is seen in absorption. How- 
ever, absorption-line surveys can measure the aver age in- 



sion line fitting, and comparison to shock models con- 
firmed that more shock-like line ratios were associated 



tegrated covering factor of the neutral gas {Vt/Aii] § |3.4.1 ) 
by observing many different systems. As seen through 
nuclear lines of sight, the integrated covering factor in- 
creases with increasing star formation rate. Systems with 
star formation rate (SFR) ^ 1 M^t^ yr~-^ show ^/Att ~ 0.1 



Chen et al. 2Q1Q[ [Bouche et al.||2Q12 ), n earby LIRGs 



SFR = 7 - 70 M0 yr-^) have VL/At: - 0.4 ( |Rupke et al. 



2005cD , and nearby ULIRGs (SFR > 70 Mc^ yr"^) have 
nyiTr - 0.8 (|Rupke et al.[2002[|2005c| |Martin|20"Q5l ). (In 

the model of a bipolar wmd, these correspond to open- 
ing angles on one side of the disk of 60°, 110°, and 160°, 
respectively.) 

The integrated covering factors of the ionized and neu- 
tral gas can also be estimated from the SRFW, using the 
radii we assume (§ |3.4.1| ). This way of estimating rt/An 
is much less certain, especially since we assume a wind 
radius. A typical value is VL/At: ~ 0.1, if we double the 
near-side covering factor to account for a far-side outflow 
(Table [1. 

For tne neutral gas, values of Q./Air determined from 
the SRFW deprojection (0.1) are much lower than those 
based on observations of many nuclear sightlines (0.8). 
The explanation we favor for this discrepancy is the fact 
that the background continuum source declines in bright- 
ness with increasing projected radius. Thus, the neu- 
tral gas becomes increasingly difficult to detect with in- 
creasing projected radius. Furthermore, these larger ra- 
dius sightlines dominate the measured l]/47r, once depro- 
jected. Thus, the integrated covering factors measured 
from our IFS observations are lower limits to the true 
values, which are probably closer to ~0.8. 

5.2. Ionized Gas Excitation 

The excitation of the ionized gas in a wind serves as a 
signpost of galaxy-scale processes. In disk galaxies in the 
local universe, gas excitation reveals wind-driven shocks 
and ionization by the wind power source. Outflowing 
warm ionized gas is typically shock-excited in cases where 
the wind is driven by a starburst alone ( Veilleux & Rupke 
2002[ [Sharp fc Bland-Hawthorii||2Q10( |Rich et al.| |2010), 
and photoionized by the AG JN radiation held in cases 
where an AGN is al so pre sent ( Veilleux et al.|20Q3 Sharp 
' fc Bland-Hawthorn|[2QTQl ). 

In local mergers, shock ionization is also prevalent 
(iMon real-Ibero et al."200 6{ |2010| jRich et al.|[2011| |Soto 
etll^ 2012; Westmoquette et al.p012). WTthm a given 
systemTforbidden-to-recombination line ratios typically 
increase with i ncreasing radius and decreasing Ha sur 



face brightn ess ( Monreal-Ibero et al.|"2010 Westmoquette 
et al. 2012). 'Fhis is suggestive of a line tlux mixing 
model, with a varying contributio n of star formation and 
shocks determining the line ratios (|Rich et al.|2011||West- 
moquette et al. 2012). 



Merger shocks coul d be driven by tidally-induced mo- 



tions or gas outflows. Monreal-Ibero et al. 
found a correlation between line ratios an( 



2006 2010) 



gas veloc- 



with higher velocity dispersion (Rich et al. 2011). Soto 
|et al.| (2012) identified two subclasses of shock- (or pos- 
sibly AGN-excited) gas - a subset with lower velocity 
dispersions (<j < 150 km s~^) which is consistent with a 
tidal origin and a subset with higher velocity dispersions 
(<j > 150 km s~^) that is probably outflowing. 

From our data, the [N II] /Ha and [O I] /Ha line ra- 
tio maps for the galaxies in our sample (e.g.. Figure 14) 
are consistent with these previously determined relation- 
ships. In individual objects, line ratios in the rotating 
gas component anticorrelate with Ha flux, correlate with 
FWHM, and increase with increasing projected galacto- 
centric radius. Within each system, small regions can 
deviate from these trends; careful inspection of the maps 
of each system show that the overall correlations do not 
tell the entire story. The exact nature of the spatial vari- 
ation of gas excitation in each source is outside the scope 
of this work, particularly since we probe a limited num- 
ber of strong emission lines. 
However, we ar e able to improve on the analysis of 



... . vv t^J. c_aj kyj.v_y u 1.1.1. a. yj±. \j \ \j 1.1. uj.J.V-y c_ajj. uiu \j 1. 

(|2012[) by studying individual systems, and 
on the analysis of Rich et aL] ([2011 ) by being able to 



ity dispersion and argued for a tidal origin of merger 
shocks. High spectral resolution, multiple profile emis- 



identify outflows. I^'igure |49| shows the ionized gas veloc- 
ity (for outflowing components) as a function of FWHM, 
color-coded by [N II] /Ha emission-line ratio. In a sim- 
ple mixing model between stellar photoionization and 
shock excitation and/or AGN photoionization, increas- 
ing the percentage of line flux from shock excitation or 
AGN photoionization to a particular spaxel will increase 
[N II] /Ha. 

We find that, within a given system, spaxels with a 
given FWHM show increased line ratios with increasing 
velocity ('^50%) with respect to systemic. Correspond- 
ingly, spaxels with a given velocity show increased line 
ratios as FWHM increases. It is clear that bulk motions 
(as traced by '^^50%) and unresolved bulk or randomized 
motions (as traced by FWHM) are both responsible for 
shock ionization of the broad emission line components in 
ULIRGs. In some systems (F10565+2448 and Mrk 273) 
these motions are due almost completely to powerful out- 
flows. In others (VV 705:NW), some of the motions 
(^50% ^ -50 km s-^ and FWHM > 400 km s"^) are 
clearly due to outflows and some (|^'5o%| ^ 50 km s~^ 
and FWHM < 400 km s~^) may be tidally induced mo- 
tions. 

In the three galaxies in our sample with lumi- 
nous AGN, photoionization by the AGN could in 
principle contribute to the observed line ratios. In 
F08572+3915:NW, the SDSS spectrum shows no evi- 
dence of NLR emission, though the AGN in this source is 
heavily obscured. In Mrk 231, there is presently no pub- 
lished spatial information on extended high-ionization 
emission lines. Finally, in Mrk 273, there is AGN- 
photoionized g as with low velocities near the SW nucleus 
(Colina et al. 1999) . However, the contribution of the 
AGJN to th e excitation of the outflowing gas in Mrk 273 
is unclear ( Colina et al.||1999| . Future IFS observations 
of these sources in high-ionization lines will help to clarify 
the origin of ionization. 

5.3. The Multiphase Wind 
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The present dataset permits the first spatiahy-resolved 
comparison between the neutral and ionized gas phases in 
merger winds. Previous spatiahy-resolved surveys have 
focused on either the ionized ga s (Heckman et al.|[l99Q 



Lehnert fc Heckman"1996^ [Veilleux et al. 2003; Sharp & 
Bland-Hawthorn 2010; Westm oquette e t al. 2012; Soto fc 
Martin 2012; Soto et al. 2012) or the neutral gas jMar^ 



tinj^OOG). Single aperture surveys have produced evi- 
dence that the velocities of the neutral and ionized gas 
phases in merger winds are not identical or correlated on 
a system-to-system basis, except perhaps in the ca se o f 
AGN (Heckm an et al.||2000|JRupke et al.||2005a|c[ ). [Soto] 
|et_ah (2012) argue that the ionized gas in ULIKG winds 
is ol comparable mass to the neutral gas. 

In Figures 50 and [5l] we show that there is a signifi- 
cant correlation between the (projected) neutral and ion- 
ized gas velocities within a given galaxy for three of the 
five galaxies in our sample with both ionized and neutral 
outflows: F10565+2448, Mrk 273, and VV 705:NW. This 
correlation is present whether '^50%, '^98% 7 FWHM is 
considered. Furthermore, in all of the systems with both 
ionized and neutral outflows, the phase with higher ve- 
locity varies from point to point and object to object. 
In 2/5 systems with both ionized and neutral outflows, 
the ionized gas is typically higher in velocity; in 2/5 the 
situation is reversed; and in 1/5 the velocities are com- 
parable, on average. 

The correlations we observe point to a physical con- 
nection between the ionized and neutral outflow phases 
within a given galaxy, even if this connection is complex. 
This complexity relates to unresolved substructure along 
the line of sight that may be different for each phase, as 
well as the fact that the tracers of these phases reflect the 
underlying gas properties in different ways (e.g., density: 
emission lines scale as and optically thin absorption 
lines as n). 

Two systems, F08572+3915:NW and Mrk 273, show 
the highest velocities in our survey, 2000 — 3000 km s~^. 
Velocities this high are seen only in the ionized gas, and 
not the neutral gas (or molecular gas; see below). These 
high velocities are also only seen in galaxies with an 
AGN, a point to which we will return in the next sec- 
tion. 

The ratio between the mass, momentum, and energy 
contained within the ionized and neutral phases varies 
from object to object (Table [t]), but in all systems 
the neutral gas dominates the mass. In two systems 
(F08572+3915:NW and Mrk 273) the ionized gas con- 
tributes significantly to the momentum and energy. The 
ionized gas, however, is clearly affected by extinction in 
these systems. Correcting for extinction would only in- 
crease the contribution of the ionized gas to a system's 
mass, momentum, and energy. It is unclear how extinc- 
tion corrections would affect the typical velocity of the 
gas. IFS observations of other Balmer lines in these sys- 
tems (in order to correct for extinction), or of ionized gas 
in the NIR and MIR, would provide better constraints on 
the ionized gas properties of the wind. 

The wind in F10565+2448 is also dusty, as shown by 
the correlation between continuum color and Na I D 



lower covering factor and lower optical depth reduce the 
absorption line equivalent width. If the gas and dust are 
closely associated, then a lower covering factor makes the 
continuum bluer by allowing more light to emerge dust- 
free, and a lower optical depth makes it bluer by simply 
reddening the continuum to a lesser degree. However, the 
underlying stellar population can also vary with position, 
and will introduce scatter and/or slope changes (within 
a given system) or systematic offsets (on a system-to- 
system basis). 

The other galaxies in our dataset do not show the 
large dynamic range in HST F435W/F814W color and 
Na I D equivalent width that are seen in F 10565 +2448. 
F17207— 0014 comes the closest, and shows a less sig- 
nificant correlation, with more scatter. The data for 
the other galaxies in our sample (except for Mrk 231, 
whose continuum structure is too strongly impacted by 
the central source for accurate colors) do not show cor- 
relations, but do lie on or near the correlation defined by 
F10565+2448. 

If the neutral phase of these winds is indeed dusty, 
then this outflowing dust is partly responsible for ob- 
scuring the visible galaxy continuum. Most of the UV 
light in ULIRG s is heavily obscured by concentration s 
of nuclear dust ( [Soifer et al.||2000| [Goldader et al.||2QQ2D ; 
however, on larger scales outnowing dust may play an im- 
portant role as a foreground screen in extinguishing the 
visible continuum. This effect is seen at high spatial res- 
olution in the kpc-scale dust filaments in F10565+2448 
(Figure 17). Future detailed analyses of the optical colors 
and other extinction probes (e.g., the Balmer decrement) 
will better illuminate this connection in other systems. 

The winds in ULIRGs a l so have a molecular compo- 
nent jFischer et al.||20 10a; "Feruglio et al."'2010! 'Sturm| 
et al |2011[|Aalto et al.|20 12; Cicone et al. 2012, Veilleux 



equivalent w idth on a spaxel-by-spaxel basis (Shih & 



Rupke 2010 Figure 52) an d the dust filaments emerg- 
ing along the minor axis (§4.2.2 and Figure 17). Both 



et al. 2013, in prep.). Except for the case of Mrk 231, 
the only information at present on the molecular phase 
in these systems comes from single-aperture Herschel ob- 
servations of OH transitions. Three galaxies in our sam- 
ple (F08572+3915:NW, Mr k 231, and F17207-0014) are 
in the 6-galaxy sample of Sturm et al.| ( |2011| ). The 
single- aperture integrated velocity profiles are not di- 
rectly comparable to our spatially-resolved data. How- 
ever, it is clear that the high neutral gas velocities in 
F08572+3915: NW and Mrk 231 ar e also observed in the 
molecular gas. Sturm et al. (2011') measure '^50%= —700 
and —600 km s~"^ and Vq8% = —1300 and —1200 km s~^, 
respectively, in these galaxies. These v^q% {vgs%) val- 
ues lie somewhere in between the average and maxi- 
mum v^Q% {vgs%) values that we measure in the neu- 
tral gas for these two galaxies, suggesting that the neu- 
tral and molecular gas are connected. However, the ion- 
ized gas in F08572+3915 reaches much higher velocities. 
In F17207-0014 (^50%= -100 and maximum velocity 
of —370 km s~^ in the molecular gas), the neutral gas 
shows significantly higher velocities than the molecular 
gas, though the highest velocity neutral gas is at large 
radius. 

The presence of molecular outflows alongside minor- 
axis neutral and ionized flows in our IFS sample suggests 
a helpful picture. On kpc scales, the dusty, molecular, 
neutral, and ionized wind emerges along the minor axis 
of a molecular disk. The molecular gas is entrained as 
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the wind emerges from this disk, and is roughly co-spatial 
with the neutral and ionized gas at small scales. 

5.4. Wind Power Sources 

Perhaps the most important outstanding issue in our 
understanding of ULIRG winds is the relative impor- 
tance of a starburst and AGN in accelerating the wind 
in a particular galaxy. This is an important question be- 
cause feedback from AGN and/or stars are frequently 
invoked in models to explain the emerg ence of QSOs 
from dusty, major mergers of disk galaxies (Sanders et al. 
1988) and the subsequen t truncation of star formation 



[Tg 

and black hole activity (Springel et al. 2005 Hopkins 



[et al.^200 5 ) . Attention has recently turned to AGIN feed- 
back in particular in this context, due to its higher ex- 
pected velocities compared to stellar feedback and its 
rapid emergence after the onset of AGN activity. 

Indirect arguments for the presence of QSO-mode feed- 
back in mergers abound, but until recently data have 
been slim. There is now evidence from multiple gas 



ample of such feedback (|Rupke et al.|2005a 


Fischer et al. 


2010a^ "l^eruglio et al.||2010| 


ISturm et al. 


2011 


IKupke 


FVeilleux 2011; Aalto et al. 


20121 Cicone et aJ 


.(|2012p. 



starburst-driven relies on the relative power of the cen- 
tral engine compared to the starburst, and on the re- 
markably high velocities observed, which are higher than 
those seen in ULIRGs with lower luminosity AGN or no 
AGN dSturm et al.|[2QlT| [Rupke fc Veilleux][2QTT ). 

Surveys of outflows m mergers and other systems based 
on nuclear spectra suggested that starbursts are suffi- 
cient to power large-scale outflows in systems where both 
are present, and found that there was no strong evi- 
dence f or the AGN dominating la rge-scale outflow dy- 
namics (Rupke et al.| |2005a |Krug et al. 2010). However, 
spatially-resolved and multiphase data, as in the case 
of Mrk 231, and the advent of precise n uclear AGN lu- 
mino sities for individual major mergers ( [Veilleux et al.| 



2009 ) have the potential to alter this conclusion. Nuclear 
Herschel spectra of nearby ULIRGs (including Mrk 231) 
show a correlation o f molecular o utflow velocity with 
AGN luminosity (Sturm et al.|2011, Veilleux et al. 2013, 
in prep.). This suggests that these (unresolved) molec- 
ular outflows are powered primarily by an AGN at the 
highest velocities. 

The current sample contains three systems whose in- 
frared luminosity is clearly dominated by star formation, 
and for which there is no clear evidence of AGN ac- 
tivity: F10565+2448, VV 705:NW, and F17207-0014. 
Our sample also has three systems that contain lumi- 
nous AGN: F08572+3915, Mrk 231, and Mrk 273. How 
do the outflows in the systems with AGN compare to 
those without? 

5.4.1. Comparison of Wind Structures 

Structurally, the winds in galaxies with and without 
AGN are not obviously different. The data are consis- 
tent with the presence of flows along the minor axis of 
a nuclear disk in almost every system, with or without 
an AGN. Furthermore, superbubbles (F08572-h3915:NW 
and VV 705 :NW) and widespread extended outflows 
(F10565+2448, Mrk 231, Mrk 273, and F17207-0014) 



are each present in nuclei with and without AGN. Ionized 
gas excitation correlates with ^50% and FWHM in both 
types of systems, suggesting shock excitation (though 
AGN photoionization may still be present; §5.2). 

A wind perpendicular to a nuclear disk is consistent 
with the structures of starburst-driven winds in the local 
universe. However, it differs from those of the narrow line 
region (NLR) outflows in nearby Seyfert galaxies, which 
have random orientations with respect to the galaxy disk 
jCrenshaw fc Kraemer"2000^ 'Crenshaw et al. 2000'; 'Fis-| 
[cEer et al. 2010b, 2011; Miiller- Sanchez et al. 2011).^ 
collimated outflow perpendicular to a galaxy disk is thus, 
at face value, more consistent with a starburst-driven 
wind than a canonical Seyfert NLR outflow. However, 
this inference is based on the study of local disk galax- 
ies, rather than mergers with massive molecular disks 
and luminous AGN. The AGN in our sample, with lumi- 
nosities 3 - 10 X 10"^^ erg s"^ | Veilleux et al.||2009| ), are 
much brighter than Seyferts. They are in fact radio-quiet 
QSOs, whose narrow line region sizes can reach severa l 
kpc or more (Bennert et al. 2002j Greene et al. 2011). 
(NLR outflows in Seyferts are typically confined to sub- 
kpc scales.) The structure of NLR outflows in radio-quiet 
QSOs is not well understood. From our data, we cannot 
determine the exact relationship between the outflowing 
ionized gas in our sample and the classical NLR, since we 
do not observe high-ionization emission lines and previ- 
ous data is inconclusive (§ |5.2| ). 

It is entirely plausible that the nuclear molecular disk 
could, on sub-kpc scales, colli mate and redirect an out- 
flow driven by the AGN (see 



Veilleux et al. 2005, and 



references therei n). Examples in the lo cal universe in- 
clude NGC 3079 ([Cecil et al.|2001[ [2002] ) and NGC 4388 
(Veille ux et al.||20'03 '). In fact, the current data point to 
exactly this occurmg in the AGN in our sample, if these 
outflows are in fact AGN-driven. NIR IFS with adaptive 
optics would be able to probe the physics of launching 
and collimation at smaller scales. 

5.4.2. Comparison of Mass, Momentum, and Energy 

We computed the mass, momentum, and energy in the 
o utflo ws in our sample using the SRFW and BSB models 
(§3.4 and Table |6|. The mass and momentum in the 
wind do not differ significantly between galaxies with and 
without AGN. However, the energy outflow rates in the 
galaxies with AGN are a factor ~4 higher than tho se with 



no AGN, primarily because of higher velocities ( §|5.4.3 ). 

The fraction of the momentum and energy in tne ion- 
ized phase (compared to the neutral phase) rises to over 
half in two o f the AGN, and is negligible in the rest of 
our sample (§[53] and Table [t]). A gain, this is due to the 
high velocities reached in the ionized gas components of 
these two systems. 

Table [8] shows the energy outflow rate in the wind com- 
pared to the mechanical luminosity of the starburst in 
each system. The mechanical lu niinosity follows f r om th e 
continuous starburst models of Leitherer et al. (1999), 
with a correction applied for a Salpeter IMl^' with lower 
mass cutoff of 0.1 Mq. In the model of starburst winds 
driven by mechanical energy from stellar winds and su- 
pernovae, there is sufficient mechanical energy from the 
starburst in each system to power the outflow. For the 
galaxies without AGN, the efficiency with which this me- 
chanical energy must be converted into the kinetic energy 
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of the outflow is low: 3% to 13%. However, for the AGNs 
in our sample, it is much higher: 100% for F08572+3915 
and Mrk 231, and >20% for Mrk 273 (with only the 
ionized phase accounted for in the latter case). A high 
thermalization efficiency is inferred for M82 (Strickland 



|& Heckman 2009), though much of the kinetic energy in 
this system may reside in the hot. X-ray emitting wind 
fluid rather than the ionized and neutral gas components. 
The fraction of the starburst mechanical energy that goes 
into the neutral and ionized phases is certainly less than 
the total thermalization efficiency, and quite possibly 
much less (Strickland & Stevens] |2000|. The starburst 
thermalization efficiencies required tor the winds in the 
AGN in our sample are thus uncomfortably high. 

If instead the wind is driven by radiation pressure from 
the hi^h luminosity in these systems (e.g.. King 2003 



Murray et al. 2005), then we can c ompare the available 



momentum injection rate (tLir/c; Murray et al. 2005| 
to the wind's momentum outflow rate. Table [s] and Pi^- 



53 show that these numbers are comparable - i.e.. 



ure 

r must be near unity for these outflows to be driven by 
radiation pressure. Given that most of the luminosity of 
these systems emerges in the far-infrared, almost every 
photon emerging from the nuclear power sources inter- 
acts at least once with surrounding interstellar dusty gas. 
Thus, radiative driving with r = 1 is a natural, though 
certainly not the only, explanation for the powerful out- 
flows in these systems. The luminosity is injected at very 
different scales for the starburst and AGN, however, since 
the MIR sizes of the AGN in our sample are much smaller 
than the system s without AGN ( Soifer et al.|2000" Iman- 
ishi et al.|[2QTT] ). 

Any radiation pressure on the wind results from both 
a starburst and any AGN. Thus, in cases where both are 
present, it is important to understand how much each is 
contributing to the outflow. Figure [53] shows that simply 
intercepting the radiation from the starburst in these sys- 
tems requires the wind to be quite optically thick (r > 4) 
in at least two cases. In the radiation pressure model, it 
is simpler to assume that the entire luminosity of the 
system, from both starburst and AGN, are powering the 
outflow. 

If an AGN drives the highest velocity flows in our sam- 
ple, then the energy outflow rate in the wind is a small 
fraction of the AGN's luminosity (0.1 — 0.4%), comfort- 
ably with in the requi rements of recent models of AGN 
feedb ack (iHopkins & Elvis 2010; Faucher-Giguere et al. 



2012|. Table ll shows the ratio ot the energy outflow rate 



in the wind to the AGN luminosity in each system in our 
sample. 

In summary, we reach several conclusions by studying 
the mass, momentum, and energy in these winds. (1) 
The outflow mass in our sample is independent of the 
galaxy's power source, within the errors. (2) The star- 
bursts in our sample are consistent with driving outflows 
via radiation pressure or mechanical energy from stel- 
lar winds and supernovae. (3) The data are consistent 
with AGN playing a significant role in driving the out- 
flow in those systems where an AGN is present (though 
the mass, momenta, and energy do not conclusively show 
that an AGN is contributing). The outflows in the AGN 
in our sample are overall more energetic, have a more 
energetic ionized phase, and require uncomfortably large 
thermalization efficiencies if they are driven solely by me- 



chanical energy from stellar winds and supernovae. In 
terms of momentum and energy, however, they are con- 
sistent with radiative driving by both the starburst and 
AGN or with recent models of AGN feedback. 

5.4.3. Comparison of Velocities 

The energetic winds in the AGN in our sample are 
largely a result of their large velocities compared to the 
starbursts. Table [3] shows the spatially- averaged and 
maximum values of '^50% and VgQ% for each gas phase 
and each galaxy, and Figure [54] shows how these values 
depend on the AGN luminosity. We show the AGN lumi- 
nosity computed from the average AGN fraction of the 
bolometric luminosity from six MIR diagnostics, as well 
as plausible uncertainty in this value using the range of 
luminosities compute d from individual MIR diagnostics 
( jVeilleux et"ar||2009[ ). We assume that the AGN lumi- 
nosities are upper limits for the galaxies in our sample 
without evidence for AGN at other wavelengths. 

Both galaxies with and without an AGN host fast out- 
flows. In systems without an AGN, we find spatially- 
averaged '^50% ('^98%) values of —250 (—450 to —650) 
km s~^ in the neutral gas and —100 (—500) km s~^ in 
the ionized gas. The highest velocities (maximum Vgg%) 
in each system (considering all spaxels and phases) range 
from —900 km s~^ to —1250 km s~^. 

In systems with an AGN, we find spatially- averaged 
^50% (^98%) values are -100 to -400 (-450 to -900) 
km s~^ in the neutral gas and —200 to —1500 (—700 to 
—2800) km s~^ in the ionized gas. The highest velocities 
(maximum '^^98%) reached in any spaxel in each system 
range from —1450 km s~^ to —3350 km s~^. 

We conclude from this analysis that the presence of 
a QSO in major mergers yields higher outflow ve- 
locities. The two systems with the highest neutral gas 
velocities (Mrk 231 and F08572+3915) contain QSOs, 
and the two systems with the highest ionized gas ve- 
locities (Mrk 273 and F08572+3915) also contain QSOs. 
Though there is clearly overlap between the velocity dis- 
tributions of the galaxies with and without AGN, only 
those with an AGN reach outflow velocities exceeding 
-1200 km s-\ and every system with an AGN exceeds 
this velocity in at least one gas phase. Furthermore, 
only the ionized gas phase reaches velocities of 
>2000 km s~^, and only in mergers hosting a 
QSO. The same trend of outflow velocity with AGN 
luminosity is seen in the molecular phase, again with 



velocities lower than 2000 km s (Sturm et al. 



2011 



Veilleux et al. 2013, in prep.), providing independent 
confirmation of this result. 

How does deprojection affect these conclusions? In 
Mrk 273 and VV 705 :NW, we can deproject in some areas 
of our data using the BSB model. Figure [55] shows how 
this deprojection affects the derived average velocities. 
In VV 705:NW, the system with the highest observed 
velocities in a non-AGN, the velocities change very little 
because the superbubble is observed face-on. In Mrk 273, 
however, the superbubble is highly inclined, and the de- 
rived velocities double. Thus, in these two sources, at 
least, deprojection only strengthens the above conclu- 
sions. 

This conclusion differs from that 



of [Rupke et al. 
( .2005aJ and Krug et al.| ( |2010J , who argued that the 
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AGN in LIRGs and ULIRGs did not accelerate the wind 
to higher velocities. However, these previous studies re- 
lied on single-aperture spectra and a single gas phase 
(neutral), as well as incomplete information about the 
AGN content of individual systems. It is clear that com- 
bining spatially-resolved, multiphase information, as in 
the current study, with more precise AGN diagnostics is 
key to fully understanding the power sources of winds in 
these galaxies. 

The maximum velocities we observe in the neutral 
phase of ULIRG winds are a lso higher than previously re- 
ported jH eckman et al |2QQQ| |Rupke et al.|2QQ2[ |2QQ5a|cj 
Martin|2D Q5, 2006^ It is clear from our spatially resolved 
data that this is because the highest velocities are often 
not found along the line of sight to the optical galaxy 
nucleus. 

The total luminosity of the system does not appear to 
determine the wind velocities in our sample (Figure 56 ) . 
The range of infrared luminosities in our sample spans 
only a factor of four, and the two highest velocity sys- 
tems are in the middle of the lu minosity range. This is 
consistent with previous results ( [Rupke et al.||2QQ5c|a ). 

The high velocities, masses, and energies we observe 
in ULIRGs with luminous AGN are comparable to the 
properties of outflows inf erred to arise on kpc scales in 



other radio-quiet QSOs (|Moe et al. 2009 Dunn et al. 
2010 Borguet et al. 201^] ^ I'he spatial constramts m 
other radio-quiet QSOs depend on detailed photoioniza- 
tion modeling and are thus not as direct as the mea- 
surements presented here. If they are correct, the simi- 
larity of these two results suggests that such QSO-driven 
flows are a general property of radio-quiet QSOs, whether 
found in mergers or otherwise. A more detailed census of 
outflows in radio-quiet QSOs of different types is clearly 
in order. 

Outflow velocities up to 1500 — 2000 km s~^ are also 
seen in what appear to be purely starburst objects: so- 
called Ly nian-break analogs, and compact star bursts at 



0.6 (^Tremonti et al. 2007 Heckman et al. 2011 



1 



Diamond-Stanic et al. 2012j though obscured AGN have 
not been ruled out in these systems). The highest outflow 
velocities from these sta rbursts apparently require very 



velocities irom tnese sta rbursts apparently require 
compact sta,r formation dOverzier et al. 2009 Hec 
et al.[|2011[|Diamond-Stanic et al.|[^IIT^] rTrLml^ 



iHeckman 
rs also 



contain compact star formation, with star formation rate 
surface densities of 100 — 1000 Mq yr~-^ (based on their 
infrared luminosi ties and molecular disk radii; [Downes 



Solomon 1998| ), comparable to Lyma n break analogs 
and the z ^ 0.6 compact starbu rsts (Overzier et al. 



2009 Diamond-Stanic et al. 



2012). However, the com- 



pact star formation in ULlKGs is heavily obscured and 
occurs amid massive reservoirs of dense gas, while the 
Lyman-break analogs and z ^ 0.6 compact starburs ts are 



UV-luminous and relatively unobscured (Heckman et al 



2005 



lata shows that an AGN is required in 



Overzier et al.||2009[ [Diamond-Stanic et al.||2 aT?) 

ufc 



TRGs 



to achieve the highest velocities. We speculate that the 
higher velocities observed in UV-luminous compact star- 
bursts compared to starburst ULIRGs result from the 
outflow not having to transfer energy and momentum 
into a dusty, dense, massive interstellar medium around 
the starburst. 



5.5. Feedback in Major Mergers 

An important outstanding question is the impact of 
these massive, high- velocity flows on their host galaxies 
and nuclear black h oles. These galaxies retai n a cen- 
tral molecular disk (Downes & Solomon 1998), as well 
as rapid star formation and active galactic nuclei, so any 
negative feedback has yet to be felt. 

In Table |8j we show that the mass outflow rate, while 
uncertain, is a significant fraction of the star formation 
rate. Thus, the outflow evacuates gas at a rate similar 
to the star formation rate as it propagates through the 
galaxy. The outflow mass will only increase if we include 



the m olecular phase (Feruglio et al.| 2010 Sturm et al. 
|2011|). Furthermore, the outflow is highly mass- loaded, 
meaning that most of its mass comes from gas it has en- 
trained while moving through the galaxy. Table [S] com- 
pares the mass outflow rate (measured) with the mass 
input from the starburst itself (from models; [Leitherer 



et al. 1999). These numbers suggest that, on average, 
about 75% of the wind mass comes from entrained ma- 
terial, and perhaps more when other gas phases are in- 
cluded. 

As we argued in § 5.1.1[ a comparison of the dy- 
namical timescales of the wind (^10 Myr) and merger 
(~100 Myr) suggests that the duty cycle of winds is 
high and that these outflows are relatively continuous. 
Similarly, the timescale for star formation to deplete 
the nuclear molecular gas reservoir in these systems 
(4 — 6 X 10^ M (7> for the four galaxies wit h CO inter- 
ferometric data; Downes & Solomon 1998) is of order 
30 — 40 Myr (the average star forniation rate in our 
sample is 143 M© yr~^). The timescale for the wind 
to deplete the nuclear gas reservoir is slightly larger: 
~100 Myr, based on the mass outflow rates we calcu- 
late here. Thus, the gas consumption, mass outflow, and 
merger dynamical timescales are all similar, and all three 
processes will be important in regulating star formation. 

The overall picture of star formation self-regulation is 
a complex one, but it is clear that gas outflows do play 
a role in governing star formation. Molecular gas is part 
of the dusty, neutral wind we detect here, and these out- 
flows will serve to end star formation in particular regions 
on a shorter timescale than would the exhaustion of the 
gas reservoir by star formation itself. This is certainly 
one means of negative feedback on star formation. 

Feedback to AGN accretion, however, is difficult to 
probe on these scales. Observations at smaller scales, 
comparison to simulations, and population studies may 
all be necessary to better constrain whether the fast AGN 
outffows in these systems are true feedback on the AGN 
itself, and on what timescales. 

6. SUMMARY 

In this study, we have presented deep integral field 
spectroscopy of the optical emission and absorption lines 
in one luminous and five ultraluminous infrared galax- 
ies at z < 0.06. Three of these galaxies have dust- 
enshrouded QSOs, and three have no apparent AGN. By 
fitting the Na I D absorption line and multiple Gaussian 
components to strong emission lines, we perform a spa- 
tial and spectral decomposition of the data. We are able 
to separate rotating from outfiowing or tidally-disturbed 
gas, and have detected neutral and ionized outfiows on 
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kpc scales in all six galaxies. This detailed data set has 
enabled us to study the multiphase structure and power 
sources of the winds in these systems. 
We reach a number of major conclusions. 

1. Where CO maps of the rotation of a nuclear disk 
are available (4 of the 6 galaxies in our sample) , the 
narrow component of the ionized gas matches the 
rotation of the molecu lar gas in those regions where 
they are coincident (§ |5.1.1 ). On larger scales, the 
narrow component of the ionized gas combines a 
continuation of the rotation of the nuclear disk with 
tidally-induced kinematics. 

2. In 5 of 6 galaxies, the dynamics and structure of at 
least one gas phase at projected galactocentric radii 
of 1 — 2 kpc are consistent with a bipolar outflow 
collimated by the nuclear molecular disk (§5.1.1). 
In two galaxies (Mrk 273 and VV 705:NW77TIiis 
bipolar outflow takes the form of a bipolar super- 
bubble that is detected on the basis of clear kine- 
matic signatures (Figures 29 and 41). On larger 
scales, the winds are less coQimate^T The preva- 
lence of minor axis flows on small scales, and their 
short timescales compared to the merger dynami- 
cal time, suggests they are a somewhat continuous 
phenomenon. 



3. Consistent with the results of |Martin| (|2006l), the 
winds in our systems reach projected galactocentric 
radii of at least 2 — 5 kpc, with the detectable size 
of the wind limited in most cases by either the size 
of our FOV or sensitivity. 



4. Consistent with the results of [Sharp _fc Bland- 



Hawthorii| (|2QlQD,|Rich et al.|(|2Qlip, and|Soto et aL 
( |2U12[ ), we nnd that the excitation of the broacT 
component of the ionized gas is due primarily to 
wind-driven shocks, since emission line ratios in 
this component increase with both velocity with re- 
spect to systemic and with line width (Figure 49). 
A contribution from AGN photoionization is also 
likely in those galaxies with a QSO. 

5. Within several systems, we flnd strong correlations 
between the velocities of the ionized and neutral 
gas phases on a spaxel-by-spaxel basis (Figures 50 
and [51]). The velocities are not often identical, how- 
ever, and it is not yet clear what determines the 
relative velocities of the two phases. 

6. We calculate that most of the mass, momentum, 
and energy in the wind is contained in the neutral 
gas phase (Table [t]). The exceptions are the galax- 
ies with the highest ionized gas velocities, which di- 
vide the momentum and energy roughly evenly be- 
tween the ionized and neutral phases. However, the 
contribution from the ionized phase may be under- 
estimated because of our inability to make extinc- 
tion corrections, and uncertainties in the neutral 
gas values remain because of uncertain ionization 
corrections. 

7. Spatially-resolved comparisons between the neutral 
phase of the wind and the optical continuum (Fig- 
ures 17 and 52 ) show that the outflows are dusty in 



at least some cases, and thus partly responsible for 
obscuring the optical continuum emission in these 
systems on kpc scales. 

In those systems without an AGN, the properties of 
the wind are consistent with being driven by radi- 
ation pressure or supernova mechanical energy. In 
those systems with a QSO, however, the required 
thermalization efliciencies appear too high to ex- 
plain the energetics on the basis of supernova me- 
chanical energy alone. Furthermore, if radiation 
pressure from a starburst alone drove the outflow 
in the QSOs, the wind optical depth would need to 
exceed unity by a factor of ^4. In the systems with 
a QSO, the data are consistent with the wind be- 
ing driven by a combination of radiation pressure 
and /or mechanical energy from the starburst and 
QSO (§ 



5.4.2). 



9. The galaxies with QSOs have the highest outflow 
velocities (Figure 54), reaching projected velocities 
of at least 1450 km s~^ in each system with a QSO, 
and up to 3350 km s-^ in one (F08572+3915:NW). 
Galaxies without AGN reach maximum projected 
velocities of ~1000 km s~^. We conclude that the 
QSO plays an important role in accelerating the 
wind in those systems in which a QSO is present. 

10. The highest gas velocities (2000 - 3000 km s"^) 
are reached in the ionized gas of two buried QSOs 
(Mrk 273 and F08572-h3915). Gas velocities this 
high are not seen in the neutral gas. 

11. When the ionized, neutral, and molecular phases 
are combined, these merger outflows are mas- 
sive enough to evacuate star-forming gas from the 
galaxy at a rate comparable to the star formation 
rate, suggesting that negative feedback is occurring 

(§11. 

Besides these scientiflc conclusions, the results of our 
study also show that single aperture nuclear spectra, or 
even long slit spectra, do not adequately capture the 
structural complexity of winds in major mergers. In par- 
ticular, IFS data is required to disentangle gas motions 
due to rotation, tidal motions, and outflow. Further- 
more, peak velocities measured in nuclear apertures, or 
along slits in a random direction, are not likely to probe 
the peak velocities in the outflow, which in our data often 
occur outside the nuclear line of sight. 

Though we are able to draw important conclusions 
from this sample, it is clear that its scope is limited 
by its size. Larger IFS samples of ULIRGs will put 
on flrmer footing, or modify, the conclusions reached 
here. Observations at higher spatial resolution (with, 
e.g., NIR IFS adaptive optics observations) will further 
illuminate the launching region and mechanism of these 
winds, particularly in those systems where the wind is 
poorly resolved (the ionized phases in F08572+3915 and 
Mrk 231) or where there are multiple, closely-spaced nu- 
clei (Mrk 273). In QSO-powered systems, studying the 
inner regions will also help to determine the physical 
mechanism by which the AGN helps to power the out- 
flow. Finally, observations of other emission lines (e.g., 
H/3 and [O III]) will enable us to better constrain the 
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connection of the dust to the ionized gas phase, the mass 
of the ionized gas phase, and the ionization mechanism 
of the emission hne gas. 
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Figure 1. (a) - (6): Two-color EST images of F08572+3915:NW at sizes \' x V and \^" x 15'^ Red is F8141^, blue is F435H^, and green 
is their average. The image orientation is indicated by the compass rose. The FOV of the GMOS observations is indicated by the box. (c) 
Logarithm of total GMOS flux integrated from 6100 to 6800 A, in arbitrary units. The image orientation is indicated by the compass rose, 
and the axis labels are in units of kpc. The near-infrared (NIR) continuum peak is circled. Isovelocity contours of the narrow ionized gas 
component are overplotted and labeled with velocity in km s"-*^. Velocities are relative to Zsys = 0.0584. The dashed line intersects the 
peaks in rotational velocity and has PA = 120° E of N. (d) Rest-frame equivalent width of Na I D, computed from line flts, in A. (e) — (/) 
Maps of the HST continuum flux ratio, log(F435W/F814W), with sizes and orientations identical to (a) and (6), respectively, (^f) Map 
of HST color within the GMOS FOV, rotated and smoothed with a Gaussian kernel of FWHM = 1" . (h) Signal-to-noise ratio of Na I D 
equivalent width. 
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Figure 2. (a) HST F160W image of F08572+3915:NW, within the CMOS FOV. A contour highlighting the nucleus is in red. The nuclear star 
cluster used to align the images is circled, as is the NIR nucleus. Slight defects resulting from interpolation over detector artifacts appear near 
the nucleus but do not affect centering, (b) HST/ ACS F814W image. Note the common alignment of the star cluster between this image and the 
F160W image, and the relative faintness of the nucleus in the optical, (c) F814W image, smoothed with a Gaussian kernel of cr = 0'.'25 to match 
ground-based seeing. Image contours are in red, and the NIR nucleus is circled, (d) GMOS 6100 — 6800 A image. Smoothed F814W contours are 
in red. Note that the smoothed F814W image is well aligned with the GMOS image. The location of the nucleus in the GMOS data is tied to the 
F160W image by common alignment with the F814W image. Axis labels are in kpc. 
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Figure 3. For F08572+3915:NW, maps of the centroid velocity {v^q%), FWHM, and most blueshifted velocity {vqq% = Vr^Q% — 2a) in 
each velocity component in the ionized gas (based on emission lines) and neutral gas (based on Na I D absorption). Each column shows 
a different velocity measure, and each row a different gas phase or velocity component. The ionized gas velocity component is labeled in 
parentheses (1 = narrow, rotating component; 2 = broad, outflowing component). Colors indicate velocity in km s"-*^. FWHM is corrected 
in quadrature for the instrumental profile. The contours, dashed lines, and open circles have the same meaning as in Figure^ 
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Figure 4. For F08572+3915:NW, maps of Ha flux in units of 1.0 x 10~^^ erg s~^ cm~2 per spaxel, or 1.1 x 10~^^ erg s~^ cm~2 
arcsec"^ (left column), [N II] 6583 A/Ho; flux ratio (central column) and [O l] 6300 A/Ho; flux ratio (right column) in each ionized gas 
velocity component. The ionized gas velocity component is labeled in parentheses (1 = narrow, rotating component; 2 = broad, outflowing 
component). The contours, dashed lines, and open circles have the same meaning as in Figure [l] 
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Figure 5. Distributions of '^^50%, FWHM, and ^^98% in each ionized gas component and the neutral gas in F08572+3915:NW. FWHM is 
corrected in quadrature for the instrumental profile. The vertical axis shows the number of spaxels in each bin, and the horizontal axis 
velocity in km s"-*^. 
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Figure 6. (a) - (6): Two-color /f^'T image of F08572+3915:SE at sizes \' x \' and 15'' x 15''. (See Figure [T] for more details.) (c) 
Logarithm of total GMOS flux integrated from 6100 to 6800 A, in arbitrary units. See Figure [llfor more details. Velocities are relative to 
^sys = 0.0585. The dashed line is along the major axis of the galaxy bulge, with PA = 140° E olf N. (d) — (e) Maps of the HST continuum 
flux ratio, log(F435W/F814W), with sizes and orientations identical to (a) and (6), respectively. (/) Map of HST color within the GMOS 
FOV, rotated and smoothed with a Gaussian kernel of FWHM = 1". 
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Figure 7. The same as Figure [3] but for F08572+3915:SE. 
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Figure 9. The same as Figure [s] but for F08572+3915:SE. 




Observed wavelength (A) 

Figure 10. Hq;+[N II] 6548, 6583 A spectra of two spaxels in the emission-hne outflow of F08572+3915:NW. The sohd black hue is the 
data, the sohd red hue the total flt, the solid blue line the continuum flt, and the dashed lines flts to individual lines. The different colors 
of the dashed lines represent different velocity components. 
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Figur e 11. The same as Figure [T] but for F10565+2448. Velocities are relative to Zgys = 0.0431. The dashed line is along the CO line of 
nodes pownes Solomon|1998l TA = 100° E of N). 
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Figure 12. The same as Figure [2] but for F10565+2448. 





Figure 13. The same as Figure [3] but for F10565+2448. 
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Figure 15. The same as Figure [s] but for F10565+2448. 
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Figure 16. Observed velocity (^'50% and ^'98%) vs. projected galactocentric radius in the outflowing ionized (broad component only) and 
neutral gas of F10565+2448. Each point represents one spaxel, and the position angle fE of N) of the spaxel with respect to the outflow 
center is represented by its color. The lines represent single radius free wind models (§ |3.4| ) at constant velocity for = 5 kpc and 10 kpc 
(the velocity of each model is given by the value at Rproj = kpc). The flat proflles ot ^'50% imply large wind radii (>5 kpc), and the 
strong asymmetry E and W of the disk requires asymmetric acceleration. The declining '^93% val ues ei ther place an upper limit on the 
wind radius (<10 kpc) or reflect the decrease in unresolved (turbulent) motions with radius. (See § |4.2.2| for more details.) 
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Figure 17. Color maps of the HST continuum flux ratio, log(F435W/F814W), in F10565+2448, with velocity contours overlaid: (a) 
ionized gas FWHM (contours: 500, 600, 700 km s"-*^); (b) neutral gas FWHM (contours: 400 km s"-*^); (c) ionized gas '^^98% (contours: 
-800, -700, -600 km s"-*^); and (d) neutral gas '^^98% (contours: -700, -600, -500 km s"-*^). The dashed line represents the position angle of 
the CO disk (Downes &: Solomon 1998), and the open circle the CMOS continuum peak. The relative color and morphology of the blue 
clumps to the left (INj are consistent with unobscured star forming regions in the near side of a galaxy disk, while the dusty, high- velocity 
gas to the right (S) is almost certainly a minor-axis outflow. The minor-axis outflow reveals dusty fllaments and unresolved motions of 
hundreds of km s"-*^, consist ent with dusty, shocked ga s. If thi s outfl ow is emerging along the minor axis of the CO disk, it is inclined ~20° 
away from the line of sight ( [Downes &: Solomon|1998] >. (See § |4.2.2| for more details.) 



36 



RuPKE & Veilleux 




Figure 18. The same as Figurem but for Mrk 231. Velocities are relative to Zsys = 0.0422. The dashed line is along the CO line of nodes 
pownes Solomon|199"8l PA E of N). 



Galactic Winds in Mergers 




Figure 19. The same as Figure [3] but for Mrk 231. 
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Figure 21. The same as Figure [s] but for Mrk 231. 
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Figure 22. The same as Figure |T] but for Mrk 273. The open circles are the locations of the NIR n uclei. Velocities are relative to 
2:.s7/.9 = 0.0373. The dashed line is along the approximate line of nodes of the gas disk around the N nucleus ( [Downes Solomon|1998l|Cole| 
let al.|1999| [Wilson et al.'|2008| PA ^ 80° E of N). 
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Figure 23. The same as Figure [2] but for Mrk 273. 
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Figure 24. The same as Figurels] but for Mrk 273. For this sy tern, the fo ur ionized gas components are: (1) rotating, (2) broad and redshifted, (3) 
broa d and blueshifted, and (4) isolated clouds with narrow profiles (§ |4.4.2| ). For component 2, ^93% is computed by adding, rather than subtracting, 
2c7 (§[33}. 
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Figure 25. The same as Figure^ but for 
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Figure 26. The same as Figure [s] but for Mrk 273. 
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Figure 27. Hq;+[N II] 6548, 6583 A spectra of four spaxels in Mrk 273: nearest the N and SW nuclei, in the near-side superbubble, and 
in the counter-propagating superbubble. The solid black line is the data, the solid red line the total fit, the solid blue line the continuum 
fit, and the dashed lines fits to individual lines. The different colors of the dashed lines represent different velocity components. Note the 
three components in each superbubble: the near- (blueshifted) and far- (redshifted) sides of the bubble, as well as the rotating component. 
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Figure 28. A color map of the number of emission-line components fit to each spaxel in Mrk 273. The axis labels are in kpc. 
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Figure 29. Position-velocity diagrams of the ionized gas in Mrk 273. Six columns are shown, with horizontal distance from the outflow center 
in kpc overlaid on each. Each ionized gas velocity component is shown (1/black solid circles — narr ow, rotating; 2/red open squares = broad, 
redshifted; 3/blue crosses = broad, blueshifted; and 4/green bowties = narrow, non- rotating; § |4.4.2| ), as well as the n = 2 bipolar superbubble 
model fit to the data (open black circles). Th e vert ical lines emerging from each point represent the spread in the velocity distribution (iblcr). 
Model parameters are given in Table [sj (See §§ |3.4.2| and |4.4.'3| for more details.) 
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Figure 30. Color maps of superbubble model fits to the emission line velocity fields in Mrk 273. The fit shown is for n = 2. The three 
columns show the data (left), model (middle), and their difference (right). In each panel, the solid line outlines the near-side bubble, and 
the dashed line the counter-bubble. Small black circles locate the two nuclei, while the red circle is the best-fit outfiow center, assuming 
symmetry along the outfiow axis. The numbers in the left and middle columns give the average velocity of that map within the model 
confines, while the nu mb ers in the right panel are the RMS velocity differences within the same area. Model parameters are given in 
Tablejs] (See §§[3X2] and [4X3] for more details.) 
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Figure 31. The same as Figure [l] but for VV 705:NW. Velocities are with respect to Zsys = 0.04035. 
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Figure 32. The same as Figure [2] but for VV 705:NW. 
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Figure 33. The same as Figure [3] but for VV 705:NW. 
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Figure 35. The same as Figure [s] but for VV 705:NW. 
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Figure 36. The same as Figure m but for VV 705:SE. Velocities are with respect to Zsys = 0.0407. The dashed hue is the estimated hue 
of nodes based on the current dataset (PA = —30°). The nuclear position is based on the /-band peak. 
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Figure 37. The same as Figure|2] but for VV 705:SE. Note that saturated or bad pixels are present near the nucleus in the F160W image, 
and that the star clusters do not appear in the F160W image; as a result, the F814W image is used to locate the nucleus. 
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Figure 40. The same as Figure [s] but for VV 705:SE. 
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Figure 41. Position- velocity diagrams of the neutral gas in VV 705:NW. Six columns are shown, with horizontal distance from the outflow 
center in kpc overlaid on each. The neutral gas ^'59% values are shown (black open circles), as well as the n = 2 bipolar superbubble model 
flt to the data (red open circles) and the narrow ionized gas component (small black points). The ver tical li nes emerging from each point 
represent the spread in the velocity distribution (zblcr). Model parameters are given in Table ID (See § [4^ for more details.) 
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Figure 42. Color maps of superbubble model fits to the absorption line velocity field in VV 705:NW. The fit shown is for n = 2. The 
three columns show the data (left), model (middle), and their difference (right). In each panel, the solid line outlines the near-side bubble, 
the dashed lines the expected counter-bubble (the counter-bubble is hidden behind the disk). A small black circle locates the nucleus and 
outfiow center. The numbers in the left and middle columns give the average velocity of that map within the model confines, while the 
number in the right panel is the RMS velocity difference within the same area. Model parameters are given in Table [5] (See §§[3X2] and 
|4.5.2| for more details.) 
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Figure 43. The same as Figure [l] but for F17207-001 4. Velocities are with respect to Zsys = 0.0430, and the dashed hue represents the 
position angle of the CO disk ^Downes Solomon|1998( ). 
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Figure 44. The same as Figure [2] but for F17207-0014. 
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Figure 45. The same as Figure [3] but for F17207-0014. 
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Figure 46. The same as Figure [4| but for F17207-0014. 
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Figure 47. The same as Figure [s] but for F17207-0014. 



400 



800 -800 -400 



Position Angle (E of N) 

180° 270° 




3 4 1 

Projected Radius (kpc) 



Figure 48. Observed velocity (v^q% and Vqs%) vs. projected galactocentric radius in the outflowing neutral gas of F17207— 0014. Each 
point represents one spaxel, and the position angle (E of N) of the spaxel with respect to the outflow center is represented by its color. The 
lines represent single radius free wind models (§13} at constant velocity for = 4 kpc and 8 kpc (the velocity for a given model is given by 



the value at R^j 



kpc). The observed increase of velocity with projected galactocentric radius is inconsistent with the SRFW model. 
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Figure 49. Ionized gas velocity ('^^50%) vs. FWHM for the broad emission line component, in the four nuclei in our sample with ionized 
gas outflows and resolved [N II] /Ho; data. The color of each point represents its [N II] /Ho; flux ratio. Within a given source, there is a clear 
trend for spaxels of higher velocity with respect to systemic and higher FWHM to have higher gas excitation. This suggests that the wind 
is being shock-excited, though there may be a contribution from AGN photoionization in those systems where an AGN is present, as well 
as a contribution from stellar photoionization at the lowest velocities (§|5.2|. 
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Figure 50. Ionized vs. neutral gas velocity, for the 5 nuclei in our sample with both ionized and neutral winds. Each point represents a 
spaxel with outflowing gas. The black points show '^^50%, and the red points ^^98%- The solid black line denotes equality. Within at leas t 
three systems (F10565+2448, Mrk 273, and VV 705:NW), there is a clear correlation between the velocities of the two gas phases (§ |5.3| >. 
It is unclear, however, what causes the deviations from equality, and why these deviations differ from object to object. 
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Neutral Gas FWHM (km/s) 

Figure 51. Ionized vs. neutral gas FWHM, for the 5 nuclei in our sample with both ionized and neutral winds. Each point represents a 
spaxel with outflowing gas. The solid black line denotes equality. Within at least three systems (F10565+2448, Mrk 273, and VV 705:NW), 
there is a clear correlation between the velocities of the two gas phases (§ |5.3| ). It is unclear, however, what causes the deviations from 
equality, and why these deviations differ from object to object. 
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Figure 52. Rest-frame equivalent width of Na I D vs. HST color. The colors are binned and aligned to correspond to GMOS spaxels. 
The dashed line is a fit to F10565+2448. The correlations for F10565+2448 and F17207-0014 suggest that the dusty, neutral wind is 
attenuating the galaxy's c ontin uum light. Though they have smaller dynamic range, the other systems fall near this correlation; these 
winds may also be dusty (§|5.3|. 
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Figure 53. Momentum outflow rate normalized to infrared luminosity and starburst luminosity vs. infrared luminosity. Large, black 
symbols represent c dp/dt/ Lipi, and small, red symbols represent c dp/dt/ Lsb- Circles (hourglasses) are results from the SRFW (BSB) 
model. The dashed lines represent required optical depths of the outflow to photons if all of Ljpi or Lsb is scattered by the wind. Note 
that powering the outflows with the starburst alone requires higher optical depths in some cases. 
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Figure 54. Velocity vs. AGN luminosity. Spatial averages and maximal values are computed from the projected velocities of outflowing 
components in all spaxels with outflowing gas. Blue (red) symbols represent neutral (ionized) gas. Two ionized gas symbols are shown for 
Mrk 273, representing the approaching and receding sides of the outflow. AGN luminosity is computed from the infrared luminosity and 
the average of six mid-infrared (MIR) diagnostics (Veilleux et al. 2009); the horizontal error bars represent the range of AGN luminosities 
determined from individual MIR diagnostics. The AGIN luminosities for systems without clear multiwavelength evidence for an AGN are 
shown as upper limits. As labeled by the arrows in the upper right panel, the dashed line divides systems with and without clear evidence 
for an AGN. It is clear that galaxies with a luminous AGN (or QSO; F08572+3914:NW, Mrk 231, and Mrk 273) reach higher outflow 
velocities (up to 1450 — 3350 km s~-^ in each object with a QSO; § |5.4.3| . 
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Figure 55. Velocity (^^50%) vs. AGN luminosity for the two systems with a bipolar superbubble (Mrk 273 and VV 705:NW), to illustrate 
projection corrections to Figure[54] Hourglasses show average deprojected velocities from the BSB model, and circles show average projected 
velocities. To enter these averages, spaxels must fall within the superbubble boundary and contain significant absorption or emission. The 
average velocities in VV 705:N W change little on deprojection, since it is a face-on bubble, while those in the highly inclined Mrk 273 
bubble roughly double (§ |5.4.3| ). 
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Figure 56. Velocity vs. total infrared luminosity (see Figure [54] for more details). There is no obvious relationship between velocity and 
infrared luminosity. The dynamic range in luminosity is only a tactor of four, and the starburst and starburst+AGN systems are mixed; 
the sources with AGN are not the brightest. 
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Table 1 

Sample 



Galaxy 


Other name 


z 




125/ J 60 


AVoriN rrac. 


or tx 


Spec. Type 


Merger Class 


iNuc. oep. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


F08572+3915 




0.0584 (NW) 
0.0585 (SE) 


12.10 


0.241 


0.72 


61 


H/C (NW) 
H (SE) 


Illb 


5.0 


F10565+2448 




0.0431 


12.02 


0.105 


0.17 


150 


C 


triple 


triple 


Mrk 231 


F12540+5708 


0.0422 


12.51 


0.287 


0.71 


162 


SI 


IVb 


0.0 


Mrk 273 


F13428+5608 


0.0373 


12.14 


0.105 


0.46 


128 


S2 


Illb 


1.0 


VV 705 


F15163+4255, 


0.04035 (NW) 


11.89 


0.157 


0.10 


120 


C (NW) 


Illb 


7.5 




Mrk 848, I Zw 107 
















0.0407 (SE) 










C (SE) 






F17207-0014 




0.0430 


12.39 


0.050 


0.11 


376 


H 


IVa 


0.0 



Note. — Col. 1: Galaxy name. Names starting with "F" are from the Infrared Astronomical Satellite (IRAS) Faint Source Catalog. Col. 2: 
Other common names. Col. 3: Redshift (geocentric), from the current study. Col. 4: 8 — 1000/im infrared luminosity, based on IRAS or Spitzer Space 
Telescope fluxes from, in order of preference: Veilleux et al. (2009); Surace et al. (2004); Sanders et al. (2003); or NED. Col. 5: 25-to-60 fim flux density 
ratio. Col. 6: "AGN fraction", or fraction of the bolometric luminosity of the galaxy due to an AGN. This value is determi ned from M IR diagnostics 
d Veilleux et al.|20Q9| ) except in the case of VV 705, for which we assume a value of 10% (§[4?5]). For ULIRGs, Lboi ~ 1-15Lir ( [Sanders Mirabel 1996] 

f!ineiD^^T'2009T. Col. 7: Star formation rate in M,^ yr~^, as determined from the AGN fraction, the infrared luminosity, and the pre script ion of 
3nnicutt (1998). Col. 8: Optical spectral type, from [Yuan et al. (2010) except for F08572+3915; we determine types for this system in § |4.1.1| H = 
ii-region-like, C = composite, Sl/2 = Seyfert 1/2. Col. 9: Merger class based on system from Veilleux et al. (2002). Classes are taken froni this paper 
for F08572+3914 and Mrk 231. For ot her sources: F 10565+244 8 class was changed fro m Ilia to triple based on archival HST images; Mrk 273 class 
was changed from IVb to Illb b ased on|Scoville et al.t (200Q) an d|Iwasawa et alT] (|2Qllb|); and VV 705 and F17207— 0014 were newly classifled. Col. 10: 
Nuclear separation in kpc, from [Scoville et^ir]'|^n0Q|) except for F10565+2448 and VV 705. The separation for VV 705 is based on archival HST data. 



Table 2 

Observations 



Galaxy 


Dates 


texp 


PA 


Dim 


(1) 


(2) 


(3) 


(4) 


(5) 


F08572+3915:NW 


10aprl4, lOmayll, 10mayl9 


3 X 1800s 


81° 


16 X 22 


F08572+3915:SE 


10apr03, 10aprl4, 10mayl9 


4 X 1800s 


81° 


16 X 24 


F10565+2448 


07febl2 


7 X 1800s 


90° 


22 X 26 


Mrk 231 


07jull9 


5 X 900s 


35° 


21 X 25 


Mrk 273 


10junl3 


4 X 1800s 


0° 


20 X 27 


VV 705:NW 


10aprl4, 10apr20, 10apr22 


5 X 1800s 


35° 


16 X 22 


VV 705:SE 


10aprl4, 10apr20, 10apr22 


4 X 1800s 


35° 


16 X 22 


F17207-0014 


10may09, lOmaylO, 10jun03 


4 X 1800s 


90° 


22 X 30 



Note. — Col. 1: Galaxy name. Col. 2: UT dates of observations (in YYmmmDD 
format). Col. 3: Number of exposures X length of each exposure, in seconds. Col. 
4: Position angle, E of N, of long axis of flnal data cube. Col. 5: Size of flnal data 
cube, in units of 0'.'3 spaxels. 



Table 3 

Outflow Velocity Statistics 



Galaxy 


Phase 


(FWHM> 


Max. FWHM 


(^50% > 


Max. v^o% 


('^98% > 


Max. vqs% 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


F08572+3915:NW 


neutral 


583 


583 


-403 


-658 


-898 


-1153 




ionized 


1522 


1522 


-1524 


-2052 


-2817 


-3345 


F10565+2448 


neutral 


293 


544 


-218 


-347 


-468 


-757 




ionized 


472 


785 


-133 


-328 


-535 


-894 


Mrk 231 


neutral 


613 


996 


-396 


-830 


-917 


-1455 




ionized 


527 


997 


-217 


-482 


-665 


-1141 


Mrk 273 


neutral 


439 


996 


-96 


-548 


-469 


-1395 




ionized (blue) 


760 


1553 


-356 


-1465 


-1002 


-2295 




ionized (red) 


737 


1553 


284 


1063 


911 


1734 


VV 705:NW 


neutral 


481 


872 


-242 


-578 


-651 


-1232 




ionized 


469 


705 


-93 


-213 


-492 


-802 


F17207-0014 


neutral 


454 


835 


-270 


-500 


-656 


-986 


Starburst ULIRGs^ 


neutral 


330 




-170 




-450 





Note. — Col. 2: Gas phase. Col. 3—4: Mean and maximum of full width at half maximum of velocity distribution, 
in km . Col. 5 — 6: Mean and maximum of velocity at center of (Gaussian) distribution. Col. 7—8: Mean and 
maximum of Vgs% = ^50% — 2cr. 



^ Computed from starburst-driven ULIRGs in |Rupke et al.1 ( |2QQ5c1 . This is a sample average, as opposed to the spatial 
averages in the current sample. 



RuPKE & Veilleux 



Table 4 

Single Radius Free Wind Parameters 



Galaxy 
(1) 


Phase 
(2) 


Radius 
(3) 


n/47r 
(4) 


F08572+3915:NW 


neutral 
ionized 


2.0 
2.0 


0.06/0.02 
0.18 


F10565+2448 


neutral 
ionized 


5.0 
5.0 


0.17/0.12 
0.11 


Mrk 231 


neutral 
ionized 


3.0 
3.0 


0.24/0.14 
0.06 


Mrk 273 


neutral 
ionized 


4.0 
4.0 


0.05/0.03 
0.16 


VV 705:NW 


neutral 
ionized 


3.0 
3.0 


0.14/0.08 
0.07 


F17207-0014 


neutral 


5.0 


0.10/0.06 



Note. — Col. 2: Gas phase. Col. 3: Fixed wind 
radius, in kpc. Col. 4: Integrated covering factor, ex- 
pressed as a fraction of A-k. The first number for the 
neutral phase does not include a correction for the fit- 
ted covering factor of t he abs orption line in each spaxel; 
the second does (see §|3.4.l|). 



Table 5 

Bipolar Superbubble Parameters 



Galaxy 


Phase 


Rz 




i 


PA 




n 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Mrk 273 


ionized 


l.OiO.l 


0.8±0.1 


75±5° 


-5± 10° 


1600 ± 500 


2^ 






l.OiO.l 


0.8±0.1 


75±5° 


-5± 10° 


2100 ± 500 


3^ 


VV 705:NW 


neutral 


2±1 


1^ 


6±6° 


35 ± 40° 


500 ± 150 


2^ 






3±2 


1^ 


4±4° 


35 ± 40° 


600 ± 150 


3^ 



Note. — Col. 2: Gas phase. Col. 3 — 4: Semi-principal axes of the ellipsoid bubble, in 
kpc. The z axis is the common axis shared by the two counterpropagating bubbles in each 
system. Col. 5: The inclination of the bubble's z axis with respect to the line of sight (for 
2 = 0°, the z axis and LOS are parallel). Col. 6: The position angle of the bubble's z axis 
(for PA = 0°, the z axis of the near side bubble is directed N). Col. 7: The (maximum) 
velocity of the bubble along the z axis, in km s~^ Col. 8: The power-law index of the 
bubble's velocity field: v = Vz(r/rz)^, where r is the distance from the center of the wind, 
and Vz and = 2Rz are the velocity and distance along the z axis. 
^ This parameter is fixed. 
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Table 6 

Mass, Momentum, and Energy 



Galaxy 


phase 


model 


M 


dM/dt 


P 


c X dp/dt 


E 


dE/dt 








Mq 


Mq/jt 


dyne s 


Lq 


erg 


erg/s 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


h 08572+3915:N W 


neutral 


SRFW 


7.88 


1.39 


48.98 


11.95 


56.69 


42.78 




ionized 


SRFW 


6.93 


0.88 


48.47 


11.84 


56.58 


43.05 




total 


SRFW 


7.92 


1.51 


49.10 


12.20 


56.94 


43.24 


h 10565+2448 


neutral 


SRFW 


9.10 


1.81 


49.80 


11.95 


57.17 


42.42 




ionized 


SRFW 


7.88 


0.16 


48.15 


10.00 


56.06 


40.93 




LOiai 


O-Ttr VV 


9 12 


1 82 




1 1 Qf^ 


^7 90 
<J ( . zu 


42 44 


Mrk 231 


neutral 


SRFW 


8.92 


2.25 


50.02 


12.82 


57.74 


43.68 




ionized 


SRFW 


6.97 


-0.10 


47.67 


10.29 


55.50 


41.17 




total 


SRFW 


8.92 


2.25 


50.02 


12.82 


57.74 


43.69 


Mrk 273 


neutral 


SRFW 


8.23 


0.90 


48.79 


11.07 


56.50 


41.87 




ionized 


SRFW 


8.07 


0.87 


48.77 


11.25 


56.76 


42.27 




total 


SRFW 


8.46 


1.19 


49.08 


11.47 


56.95 


42.41 




ionized 


BSB 


7.27 


1.18 


48.35 


11.66 


56.36 


42.77 


VV 705:NW 


neutral 


SRFW 


8.53 


1.59 


49.36 


11.89 


56.97 


42.61 




ionized 


SRFW 


7.66 


0.28 


48.05 


10.13 


55.92 


41.09 




total 


SRFW 


8.59 


1.61 


49.38 


11.90 


57.01 


42.62 




neutral 


BSB 


8.22 


1.33 


48.86 


11.36 


56.55 


42.16 


F17207-0014 


neutral 


SRFW 


8.63 


1.53 


49.52 


11.87 


57.05 


42.52 



Note. — Col. (2): Gas phase. Col. (3): Outflow model. SRFW = Single radius, mass con- 
serving free wind. BSB = Superbubble. Col. (4) — (9): Logarithms of the mass, mass outflow rate , 
momentum, momentum outflow rate, energy, and energy outflow rate contained in the wind (§ |3.4| ). 
The ionized gas values depend on electron density as n~^, and we assume rie = 10 cm""^. The 
neutral gas values depends on ionization fraction y as {1 — y)~^ , and we assume y = 0.9. Note that 
correcting the ionized gas for reddening will increase the values in this phase, as well as the total. 



Table 7 

Percent of Mass, Momentum, and Energy in Neutral Phase 



Galaxy 


model 


M 


dM/dt 


P 


dp/dt 


E 


dE/dt 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


F08572+3915:NW 


SRFW 


89 


76 


76 


56 


56 


34 


F10565+2448 


SRFW 


94 


97 


97 


98 


92 


96 


Mrk 231 


SRFW 


98 


99 


99 


99 


99 


99 


Mrk 273 


SRFW 


59 


51 


51 


39 


35 


28 


VV 705:NW 


SRFW 


88 


95 


95 


98 


91 


97 


Note. Col. 2: Outflow model. SRFW = Single 


radius, mass conserving 



free wind. BSB = Superbubble. Col. 3—8: Percent of the mass, mass outflow 
rate, momentum, momentum outflow rate, energy, and energy outflow rate in 
the neutral phase of the wind (compared to the total in the neutral and ionized 
phases). Note that correcting the ionized gas for reddening will decrease these 
values. 
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Table 8 

Power Sources and Feedback 



Galaxy 


phase 


model 


c X dp/dt/ 


dE/dt/ 


dE/dt/ 


dM/dt/ 


dM/dt/ 










dE/dtsB 


Lagn 


SFR 


dM/dt SB 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 




neutral 


ortr vv 


-0 21 








U. (JO 




ionized 


SRFW 


-0.32 


-0.17 


>-2.61 


-0.90 


0.08 






SRFW 


0.04 


0.01 


>-2.42 


-0.28 


0.70 




neutral 


o-TLr vv 


n 1 

-U. lO 


1 1 Q 

- 1 . It) 




-u.ou 


n fi9 
u.uz 




ionized 


o-TLr vv 


-z.uo 


-Z.Do 


A 09 


9 ni 

-Z.Ul 


1 n*^ 

-l.UO 




total 


SRFW 


-0.13 


-1.18 


-2.52 


-0.35 


0.63 


Mrk 231 


neutral 


SRFW 


0.25 


0.03 


>-2.38 


0.04 


1.02 




ionized 


SRFW 


-2.28 


-2.48 


>-4.90 


-2.31 


-1.33 




total 


SRFW 


0.25 


0.04 


>-2.38 


0.05 


1.03 


Mrk 273 


neutral 


SRFW 


-1.13 


-1.68 


>-3.64 


-1.21 


-0.23 




ionized 


SRFW 


-0.95 


-1.28 


>-3.24 


-1.23 


-0.25 




total 


SRFW 


-0.73 


-1.14 


>-3.09 


-0.92 


0.06 




ionized 


BSB 


-0.54 


-0.78 


>-2.74 


-0.93 


0.05 


VV 705:NW 


neutral 


SRFW 


-0.06 


-0.91 


-1.99 


-0.49 


0.49 




ionized 


SRFW 


-1.82 


-2.44 


-3.51 


-1.80 


-0.82 




total 


SRFW 


-0.05 


-0.90 


-1.97 


-0.47 


0.51 




neutral 


BSB 


-0.59 


-1.36 


-2.44 


-0.75 


0.23 


F17207-0014 


neutral 


SRFW 


-0.58 


-1.49 


-2.61 


-1.05 


-0.07 



Note. — Col. 2: Gas phase. Col. 3: Outflow model. SRFW = Singe radius, mass conserving 
free wind. BSB = Superbubble. Col. 4: Logarithm of the momentum outflow rate divided by the 
momentum input rate if all of the galaxy's infrared luminosity accelerates the wind and each photon 
is intercepted once. Col. 5: Logarithm of the ener gy outflow rate divided by the mechanical energy 
production rate from a continuous starburst ftLeitherer et al.|1999| . Col. 6: Logarithm of the energy 
outflow rate divided by the AGN luminosity. Col. 7: Logarithm of mass outflow rate normalized to the 
star formation rate. Col. 8: Logarithm of t he mass outflow rate d ivided by the hot gas mass production 
rate expected from a continuous starburst ( [Leitherer et aL |1999|). 



